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summarizes  the  status  of  the  21 -channel  DARPA/RADC,  visible- 
wavelength,  multidither  COAT  experimental  model  and  associated 
hardware  and  its  use  in  the  experimental  measurements  on  this  con¬ 
tract.  Computer  simulation  studies  of  propagation  distortion  com¬ 
pensation  and  complex  target  (speckle)  effects  are  also  described. 
Recommendations  are  made  for  future  work. 

Experimental  observations  with  the  21 -channel  COAT  system  show 
that  blooming  distortions  occurring  in  the  first  30%  of  the  focused 
propagation  path  can  be  compensated,  leading  to  roughly  a  factor  of 
1.  5  increase  in  peak  focused-beam  irradiance.  Depending  on  the 
experimental  geometry,  however,  correction  factors  from  1.0  to  4.0 
have  been  observed.  Both  experimental  and  computer  simulation 
results  have  shown  that  turbulence  compensation  performance  is  not 
degraded  by  the  presence  of  thermal  blooming,  even  though  little 
blooming  compensation  may  occur.  Experimentally,  the  addition  of 
tracking  and  focus  controls  to  18-channel,  planar-array,  COAT  phase 
control  had  no  effect  for  blooming  distortions,  but  produced  some 
correction  for  an  artificially-generated  turbulence  alone.  The 
18-channel  phase  controls  could  remove  essentially  100%  of  the  turbu¬ 
lence  distortions  in  most  cases  where  there  was  no  significant  wave- 
front  tilt  error. 

Computer  simulation  of  phase-conjugate  (return- wave)  COAT  systems 
has  shown  almost  no  blooming  compensation.  In  addition,  there  is 
some  evidence  that  a  return-wave  COAT  system  may  reduce  the 
focal-plane  irradiance  if  the  blooming  distortions  become  too  strong. 

Experimental  measurements  of  COAT  operation  with  equal  transmitter 
and  receiver  apertures  and  semidiffuse,  extended- glint  scotchlight 
surfaces  have  failed  to  produce  any  degradation  in  the  system  converg¬ 
ence  level,  even  though  significant  spurious  modulations  within  the 
dither  band  were  observed  in  the  COAT  receiver.  Experimental  and 
computer  simulation  results  have  shown,  however,  that  it  is  possible 
to  reduce  the  COAT  system  convergence  level  if  spurious,  multiplica¬ 
tive  receiver  noise  signals  are  strong  enough.  Such  signals  can  be 
produced  in  a  COAT  receiver  when  the  speckle  pattern  return  (the 
target  signature)  moves  by  the  COAT  receiver.  A  simple  analytical 
model  has  been  developed  that  can  accurately  predict  the  convergence 
level  of  a  multidither  COAT  system  in  the  presence  of  spurious 
receiver  modulation  noise.  These  studies,  as  well  as  others,  indicate 
that  there  are  potential  problems  with  speckle-induced  noise  in  a  COAT 
receiver.  Additional  studies  will  be  required,  however,  to  show  what 
classes  of  targets  and  scenarios  can  produce  spurious  receiver  signals 
that  can  significantly  degrade  the  performance  of  a  multidither  COAT 
system. 
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PREFACE 


This  report  was  prepared  by  Hughes  Research  Laboratories, 

Malibu,  California  under  Contract  F30602-75-C-0001 .  It  describes  work 
performed  from  1  April  1975  to  1  June  1975  and  summarizes  the  principal 
results  of  the  entire  contract  during  the  period  1  July  1974  to  1  July  1975. 
The  principal  investigator  and  principal  scientist  is  Dr.  James  E.  Pearson. 
The  project  is  part  of  the  adaptive  optics  program  in  the  Opto- Electronics 
Department,  managed  by  Dr.  Viktor  Evtuhov,  at  the  Hughes  Research 
Laboratories. 
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GLOSSARY  OF  ABBREVIATIONS  AND  ACRONYMS 


1. 

COAT: 

Coherent  Optical  Adaptive  Techniques.  Any  one  of  a 
number  of  adaptive  optical  techniques.  Also  used  as  a 
descriptive  modifier  as  in  "the  COAT  system." 

2. 

BAR: 

Beam  Active  Track.  Name  given  to  the  multidither 
tracking  and  focus  controls  built  and  studied  on  this 
contract. 

3. 

HRL: 

Hughes  Research  Laboratories,  Malibu,  California 

4. 

GRC: 

General  Research  Corporation,  McLean,  Virginia 

5. 

RADC: 

Rome  Air  Development  Center,  Rome,  New  York 

6. 

NSWC: 

Naval  Surface  Weapons  Center,  Silver  Spring,  Maryland 

7. 

DARPA: 

Defense  Advanced  Research  Projects  Agency, 
Washington,  D.C. 
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SUMMARY 


This  research  program  on  coherent  optical  adaptive  techniques 
(COAT)  has  three  primary  objectives:  (1)  to  experimentally  study  the 
ability  of  a  multidither  COAT  system  to  compensate  for  thermal  blooming 
and  turbulence,  particularly  when  it  operates  against  complex,  multiple- 
glint  targets;  (2)  to  use  computer  simulations  to  cros s-check  the  experi¬ 
mental  results  and  to  investigate  areas  beyond  the  capabilities  of  the 
experimental  hardware;  (3)  to  perform  preliminary  studies  of  the  effects 
of  speckle  noise  modulation  on  COAT  system  convergence  and  performance. 
This  report  is  the  final  contract  report  and  reviews  the  experimental, 
analytical,  and  computer  simulation  work  performed  during  the  12-month 
contract  period  from  1  July  1974  to  1  July  1975. 

The  experimental  studies  employed  the  18-channel  multidither  COAT 
planar  array,  which  operates  at  0.488  pm  and  was  constructed  on  a  pre¬ 
vious  program.  This  system  and  other  related  hardware  is  discussed  in 
Section  II  of  the  report.  Thermal  blooming  compensation  measurements 
were  made  using  this  COAT  system  with  four  different  absorption  cells: 
two  liquid  and  two  gas.  The  two  gas  cells,  each  50  cm  long,  were  con¬ 
structed  during  this  contract.  The  liquid  cells  were  part  of  the  Hughes 

program.  Each  cell  was  used  in  a  slightly  different  beam  geometry 
as  discussed  in  Section  III,  but  each  geometry  can  be  scaled  to  an  interest¬ 
ing  10.  6  pm  scenario.  The  ratio  of  the  peak  target  irradiance  with  COAT 

correction  to  that  without  was  observed  to  be  1. 0,  1.3,  1.5,  or  4.  0  depend¬ 

ing  on  the  experimental  geometry.  No  differences  were  observed  between 
blooming  obtained  with  liquid  or  gaseous  media.  The  greatest  correction 
is  observed  when  the  blooming  distortions  are  in  the  first  30%  of  the 
focused  propagation  path;  the  least  correction  occurs  when  the  blooming 
is  confined  to  the  last  50  to  60%  of  path.  In  all  cases,  the  COAT  system 
convergence  level  was  stable,  with  no  evidence  of  oscillation  or  instability 
as  the  blooming  distortions  were  increased  well  beyond  that  present  with 
the  optimum  transmitter  power. 
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Multidither  adaptive  tracking  and  focus  controls  were  designed  and 
constructed  to  be  used  with  the  18-channel  COAT  system.  The  design 
goals  of  a  100  Hz  tracking  system  and  a  10  Hz  focusing  system  were  met  by 
the  hardware.  In  the  thermal  blooming  measurements,  these  controls  pro¬ 
duced  no  further  increase  in  peak  target  irradiance  over  that  provided  by 
the  18-channel  planar  array;  in  one  turbulence-only  experiment,  however, 
some  correction  was  observed. 

We  have  experimentally  studied  the  effects  of  target  motion- induced 
speckle  modulation  on  the  18-channel  COAT  system.  A  scaled  laboratory 
experiment  using  a  rotating  scotchlite  target  produced  a  wideband  speckle 
noise  spectrum  over  the  dither  frequency  band  which  did  not  observably 
interfere  with  the  convergence  properties  of  the  system.  In  another  series 
of  experiments,  an  acousto-optic  modulator  with  64%  modulation  efficiency 
was  placed  in  the  DARPA/RADC  COAT  optical  system  and  swept  over  a 
variable  frequency  bandwidth  in  a  time  comparable  to  the  servo  loop  con¬ 
vergence  time.  The  speckle  modulation  artificially  generated  by  this  tech¬ 
nique  produced  interference  signals  over  the  dither  band  greater  than  the 
normal  dither  signals,  and  produced  converged  power  degradation  of  35%. 

Computer  simulation  modeling  of  COAT  system  performance  in  the 
.presence  of  speckle  noise  modulation  has  been  used  to  obtain  the  average 
convergence  level  of  an  18- channel  multidither  COAT  system  in  the  presence 
of  this  noise.  An  oscillation  of  the  system  about  the  average  convergence 
level  has  been  observed.  A  statistical  analysis  of  the  COAT  system  and  of 
the  speckle  noise  has  been  developed  for  use  in  explaining  the  simulation 
results.  Good  agreement  has  been  observed  between  the  analytical  model 
and  the  computer  simulation  data.  The  experimental  and  theoretical 
results  obtained  thus  far  indicate  that  the  system  performance  can  be 
significantly  degraded  if  the  modulation  is  large  enough,  but  the  systematic 
evaluation  of  which  targets,  target  ranges,  and  target  motions,  if  any,  can 
actually  produce  such  large  modulations  are  being  pursued  on  other  pro¬ 
grams.  To  date,  however,  our  experimental  results  with  the  18-channel 
COAT  system  (equal  transmitter  and  receiver  apertures)  and  realistic 
targets  indicate  that  speckle  modulations  will  not  seriously  degrade  COAT 
system  performance.  This  is  a  significant  conclusion  since  any  operational 
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COAT  system  will  deal  with  s emidiffuse,  multiple  glint  targets  that  do 
produce  speckle-like  target  signatures. 

Turbulence  compensation  in  the  presence  of  thermal  blooming  has 
been  studied  by  experiment  and  by  computer  simulation.  We  have  observed 
that  turbulence  compensation  is  not  significantly  degraded  by  the  existence 
of  thermal  blooming,  even  though  the  thermal  blooming  distortions  may  not 
be  removed  by  the  COAT  system.  This  result  is  particularly  significant  for 
laser  systems  operating  at  3. 8  pm  or  at  10.  6  pm  with  larger  (>l  m)  aper¬ 
tures  where  turbulence  distortions  dominate  the  beam  degradation.  Just  as 
blooming  and  turbulence  can  be  treated  separately  and  later  combined  to 
determine  the  total  beam  distortion,  so  can  the  COAT  compensation  for  each 
type  of  distortion  be  treated  separately.  In  the  experimental  measurements, 
identical  compensation  performance  was  observed  with  single  glint  targets 
and  with  semidiffuse,  rough- surface,  rotating  spherical  targets. 

Our  computer  simulation  studies  of  COAT  compensation  for  thermal 
blooming  have  been  restricted  to  return-wave  (e.  g.  ,  phase-conjugate)  COAT 
systems  because  of  the  large  computational  cost  when  the  multidither  servo 
code  is  combined  with  the  nonlinear  propagation  code.  This  point  is  dis¬ 
cussed  further  in  Section  IV.  The  operation  of  return-wave  systems  is  sig¬ 
nificantly  different  with  nonlinear  distortions  than  that  of  a  multidither, 
outgoing-wave  system  (like  the  18-channel  experimental  system).  We  have 
observed  very  little  blooming  compensation  in  the  simulations  and  have 
found  that  return-wave  systems  can  actually  reduce  the  target  irradiance 
when  the  transmitted  power  exceeds  the  optimum  power  (that  power  which 
produces  maximum  target  irradiance).  We  suspect  that  this  type  of 
behavior  may  be  characteristic  of  any  return-wave  system,  including  those 
that  employ  compensated  imaging  concepts  (imaging-COAT  systems). 

Based  on  previous  data*  and  on  the  investigations  summarized  in 
this  report,  we  conclude  that  a  multidither,  outgoing- wave  COAT  system 
can  remove  70  to  100%  of  turbulence-induced  beam  distortions,  100%  of 
fixed  beam  distortions  (optics,  etc. ),  can  provide  about  a  factor  of  1 . 5  to 
2.0  improvement  in  peak  beam  irradiance  in  the  presence  of  thermal 


*J.  E.  Pearson,  "Coherent  Optical  Adaptive  Techniques,"  Contract  F30602- 
73-C-0248m  Final  Technical  Report,  RADC-TR-85-46,  Jan.  1975. 
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blooming,  and  can  accomplish  all  of  these  functions  simultaneously,  limited 
only  by  the  number  of  correction  channels  and  by  the  system  speed  of 
response.  The  compensation  can  be  achieved  with  realistic  targets  having 
multiple  glints,  or  semidiffuse  surfaces.  Propagation  of  high  energy  laser 
beams  through  the  atmosphere  to  produce  large  focal-plane  irradiances 
should  thus  not  be  limited  by  the  atmosphere  if  suitable  COAT  system 
designs  are  used  and  if  severe  thermal  blooming  distortions  can  be  avoided. 
Additional  work  is  indicated  in  several  areas:  (1)  speckle  modulation  effects 
to  determine  what  targets  and  scenarios  can  produce  speckle  modulations 
that  are  strong  enough  to  degrade  multidither  COAT  performance  and  to 
design  COAT  systems  that  are  insensitive  to  such  modulations;  (2)  thermal 
blooming  compensation  —  to  investigate  techniques  that  can  minimize  thermal 
blooming  distortions;  (3)  computer  simulation  —  to  develop  codes  that  can 
simulate  thermal  blooming  compensation  with  multidither  COAT  systems 
and  realistic,  multiple- glint  targets;  (4)  COAT  operation  at  infrared  wave¬ 
lengths  —  to  demonstrate  the  hardware  technology  required  for  COAT  sys¬ 
tems  at  longer  wavelengths  and  higher  powers  and  to  investigate  target 
effects  (speckle,  glint  destruction)  with  realistic  target  surfaces  and 
geometries;  (5)  new  COAT  concepts  —  to  investigate  experimentally  new 
servo  and  system  concepts  that  offer  improved  compensation  performance, 
greater  versatility,  greater  simplicity  in  hardware  design  and  implementa¬ 
tion  and  more  reliable  operation  under  a  wider  variety  of  operational 
conditions. 
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I. 


INTRODUCTION 


A.  Program  Objectives 

There  are  two  primary  objectives  of  this  program.  The  first 
objective  is  to  determine  the  performance  limits  of  multidither  coherent 
optical  adaptive  techniques  (COAT)  by  performing  scaled  laboratory  experi¬ 
ments  which  are  designed  to  produce  quantitative  data  on  the  nature  of 
thermal  blooming  and  turbulence  and  on  the  ability  of  multidither  COAT  to 
correct  for  these  distortions,  particularly  when  the  target  has  many  moving, 
time-varying  glints.  The  second  objective  is  to  use  computer  simulation  as 
a  cross-check  on  the  experiments  and  as  an  analytical  tool  to  extend  the 
understanding  of  propagation  distortions  and  of  COAT  systems  beyond  the 
capabilities  of  the  experimental  hardware. 

B.  Research  Program  Plan 

The  performance  of  multidither  COAT  has  already  been  proved  with 
an  18- element  visible  wavelength  system  developed  and  tested  on  DARPA/ 
RADC  Contract  F30602-73-C-0248  which  was  concluded  in  July  1974.  This 
same  system  is  used  for  the  laboratory  experiments  in  this  program.  The 
computer  simulation  codes  for  atmospheric  turbulence,  thermal  blooming, 
and  the  COAT  system  were  developed  on  DARPA/NOL  Contract  N60921- 
74-C-0249.  The  previous  DARPA / RADC  COAT  contract  also  supplied  a 
design  for  a  gas  absorption  cell  which  has  been  used  to  simulate  convection- 
dominated  thermal  blooming  in  the  atmosphere. 

The  research  program  for  this  contract,  illustrated  in  Fig.  1,  runs 
from  1  July  1974  through  30  June  1975.  The  required  oral  presentation  was 
made  on  3  December  1974  as  part  of  the  DARPA/NOL  Adaptive  Optics 
Symposium  held  at  Lincoln  Laboratories.1  A  contract  amendment  negotiated 
during  the  second  contract  quarter  provided  for  the  addition  of  auto-tracking 
and  autofocus  controls  to  both  the  DARPA/RADC  COAT  hardware  and  the 
Naval  Ordnance  Laboratories  computer  code.  A  second  amendment,  nego¬ 
tiated  during  the  third  quarter,  eliminated  a  high-power  design  guidelines 
task  so  that  time  and  funds  could  be  directed  toward  experiments  on  COAT/ 
target- signature  interactive  effects.  The  schedule  in  Fig.  1  reflects  these 
changes. 
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Fig.  1.  Revised  research  program  plan  showing  COAT  tasks  and 
scheduling. 


20 


c. 


Organization  of  This  Report 


This  report  serves  two  functions:  (1)  it  presents  the  accomplishments 
of  the  fourth  contract  quarter  (1  April  1975  to  1  July  1975);  (2)  it  summarizes 
the  principal  accomplishments  of  the  entire  contract.  The  report  is  thus  both 
the  fourth  quarterly  technical  report  and  the  final  contract  report. 

Section  II  reviews  the  status  of  the  experimental  hardware  that  has 
come  to  be  known  as  the  DA.RPA/RADC  COAT  system.  Items  summarized 
include  the  18-channel  planar  COAT  array,  the  3-channel  Beam  Active 
Track  (BAT)  system,  the  associated  optical  systems  used  in  experimental 
laboratory  studies,  and  the  laboratory  distortion  generators  used  to  produce 
thermal  blooming  and  turbulence. 

Section  III  summarizes  the  experimental  results  obtained  on  turbu¬ 
lence  and  thermal  blooming  compensation.  Preliminary  investigations  on 
COAT  system  operation  with  complex,  extended -glint  targets  are  presented. 
The  section  briefly  discusses  the  effects  of  spurious  receiver  modulations 
produced  by  moving  speckle-patterns  (the  target  signature). 

Section  IV  summarizes  the  results  of  the  analytical  efforts  directed 
toward  investigating  turbulence  and  thermal  blooming  compensation  as  well 
as  complex  target  effects.  The  data  presented  were  obtained  primarily  by 
computer  simulation  of  the  propagation  distortions  and  the  COAT  system 
(either  phase  conjugate  or  multidither). 

Section  V  summarizes  the  principal  results  of  the  contract  and  the 
conclusions  that  can  be  drawn  from  this  work.  The  implications  for  addi¬ 
tional  research  that  should  follow  this  work  are  presented. 
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II. 


COAT  SYSTEM  HARDWARE 


A .  Planar  Array  Phase  Controls 

The  original  DARPA/RADC  COAT  system  was  developed  on 

Contract  F30602-73-C-0248.  A  detailed  description  of  the  system  is  pre- 

1-5 

sented  in  the  five  contract  reports  for  that  program.  A  block  diagram 
of  the  complete  system  is  shown  in  Fig.  2  and  Fig.  3  is  a  photograph  of  the 
signal  conditioning  and  servo  electronics  for  the  18-phase  control  channels. 
The  features  and  performance  of  the  system  are  summarized  in  Tables  1 
and  2.  As  can  be  seen  from  Table  2,  all  of  the  fundamental  multidither 
COAT  concepts  have  been  successfully  demonstrated  with  this  system.  The 
remaining  issues  are  no  longer  concerned  with  whether  or  not  the  system 
will  work,  but  rather  with  how  well  it  will  work  under  various  conditions. 
The  data  obtained  on  this  contract  answer  the  latter  question  for  some 
cases  and  succeeding  programs  now  in  progress  are  designed  to  further 
quantify  the  degree  of  successful  COAT  compensation  for  turbulence  and 
blooming  with  realistic  targets. 


Table  1.  DARPA/RADC  Planar  Array  COAT  System  Features 


1.8  W,  0.488  (Jim,  argon  laser 

Planar  output  arrays  with  any  configuration  containing  up  to  18  phase- 
controlled  elements  and  one  reference  element. 

PZT  cylinders  and  bimorphs  used  as  phase  dither  and  phase  correc¬ 
tion  drivers  (’’piston11  elements).  Up  to  60°  peak  dither  modulation 
and  over  ±2  X  correction  available  at  X  =  0.488  |j,m. 

Dither  frequencies  adjustable  from  8  to  32  kHz;  1  channel/frequency 
or  sine/cosine  operation. 

Chopper-stabilized  AGC  with  46  dB  dynamic  range  and  1  msec 
response  time. 

Sample-and-hold,  slew-offset  capability  —  electronic  and  mechanical. 

Two  photomultiplier  detectors  available:  one  for  ’’local  loop”  opera¬ 
tion  (no  path  distortion),  one  for  "target  loop”  operation  (path  contains 
turbulence  and  blooming). 
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(19th  ELEMENT  IN 
CENTER  IS  NOT 
SERVO-CONTROLLED) 


Fig.  2.  COAT  system  block  diagram  showing  additional  multidither  servo 
loops  for  tracking  and  focus  control. 
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Fig.  3. 


RADC/COAT  multidither  electronics:  18-phase- 
control  servo  channels  are  shown;  the  2  tracking 
(x  and  y)  channels  and  1  focus  channel  are  shown  in 
Fig.  4. 
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Table  2.  DARPA/RADC  Planar  Array  Performance 


Convergence  time  of  1.  5  ms  or  less  (10%  to  90%). 

Good  agreement  between  computer  simulation  of  servo  and  experi¬ 
mental  results. 

Single  glint  tracking  at  rates  up  to  30  mrad/sec. 

L  3  x  diffraction-limited  performance  (no  propagation  distortions). 
Convergence  on  strongest  glint  in  a  moving,  multiple  glint  target. 
Demonstrated  turbulence  compensation. 

Demonstrated  thermal  blooming  compensation. 

Demonstrated  turbulence  compensation  in  the  presence  of  blooming. 

Demonstrated  compensation  performance  with  M extended- glint11 
targets. 


B.  Multidither  Tracking  and  Focus  Controls 

Amendment  No.  1  to  the  original  contract  provided  supplemental 
funds  to  add  three  additional  multidither  control  channels  to  the  18-channel 
planar  array  controls:  2  tracking  controls  and  1  spherical-focus  control. 
Although  we  usually  will  refer  to  the  entire  COAT  system  as  the  n21  -channel 
DARPA/RADC  COAT  system,  M  we  will  have  occasion  to  refer  to  just  the 
focus  and/or  tracking  controls.  In  this  report,  we  will  use  the  term 
n18-channeln  to  refer  to  the  18-element  planar  array  system  which  pro¬ 
duces  piston-type  phase  control  on  18  separate  transmitted  beam  segments. 
The  acronym  nBAT,T  (for  Beam  Active  Tracking)  will  be  applied  to  the 
focus  and/or  tracking  controls  alone.  When  speaking  of  combinations,  we 
will  use  the  terms  M20-channelM  (COAT  plus  tracking),  or  n21-channelM 
(COAT  plus  focus  plus  tracking). 

The  tracking  and  focus  controls  are  shown  in  Fig.  4  and  Fig.  5 
is  a  schematic  of  the  focus  actuator.  The  servo  control  channels  were 
built  by  dividing  the  servo  parameters  of  the  18-channel  system  by  10  for 
the  tracking  channels  and  by  100  for  the  focus  control.  This  procedure 
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Fig.  5. 

Piezo- electrically- driven,  variable -radius 
spherical  mirror  used  for  autofocus  control. 
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was  expected  to  produce  a  tracking  servo  with  a  1  kHz  dither  and  a  10  msec 
convergence  time  (100  Hz  response)  and  a  focus  servo  with  a  100  Hz  dither 
and  a  100  msec  convergence  time  (10  Hz  response).  The  oscilloscope  traces 
of  target  irradiance  versus  time  shown  in  Fig.  6  demonstrate  the  achieve¬ 
ment  of  these  design  goals. 

Both  controls  could  be  made  to  operate  faster  since  the  lowest 
resonance  frequency  of  the  focus  actuator  is  3.4  kHz  and  the  tracking 
galvanometer-plus -mirror  resonance  is  at  1.7  kHz.  Although  we  chose  to 
operate  the  two  tracking  controls  with  two  separate  dither  frequencies 
(750  Hz  and  1  kHz),  sine-cosine  operation  will  work  equally  well  (producing 
a  conical-scan  tracker). 

The  tracking  performance  of  the  controls  is  shown  in  Fig.  7.  An 

interesting  behavior  is  observed:  the  controls  do  not  automatically  steer 

the  beam  boresight  axis  onto  the  glint  when  the  phase  and  BAT  controls  are 

used  together.  The  stable  lock  position  depends  on  the  initial  conditions 

as  shown  in  the  figures.  This  behavior  is  not  entirely  reproducible  and  is 

a  consequence  of  the  sidelobe  structure  of  the  segmented  transmitting 

aperture  and  of  the  choice  of  dither  frequencies.  (The  behavior  would  not 

exist  if  the  tracking  servo  were  faster  than  the  phase  control  servo.)  This 

explanation  is  discussed  in  greater  detail  in  the  third  contract  quarterly 

.  6 

report. 

Although  the  focus  control  servo  has  performed  as  predicted  by  the 
design,  the  actuator  shown  in  Fig.  5  has  one  problem:  it  introduces  a 
small,  but  not  insignificant  amount  of  aberration  into  the  optical  beam. 

When  the  actuator  is  used  with  the  18-channel  COAT  system,  the  distortion 
is  no  problem  since  the  COAT  system  can  completely  remove  it.  The  dis¬ 
tortion  does  not  allow  the  actuator  to  be  used  by  itself,  however,  without 
introducing  undesired  beam  distortions.  Consequently,  all  of  the  experi¬ 
ments  discussed  in  Section  III  that  involve  the  focus  control  also  use  the 
18-channel  COAT  system. 

The  BAT  tracking  and  focus  controls  have  also  been  implemented 
into  a  multidither  computer  simulation  that  uses  a  deformable  mirror  as  a 
dither  and  corrector  element.  The  nonboresighting  effects  seen  in  the 
experiment  were  thus  not  observed.  We  also  decided  to  implement  a 
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Measured  convergence  time  of  BAT 
controls,  (a)  Tracking  controls 
(b)  Focus  control. 
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Fig.  7. 

Tracking  performance  of  COAT  controls,  (a)  Tracking  control 
off.  The  formed  array  pattern  is  also  shown  for  reference. 

(b)  Tracking  control  performance.  Curves  with  tracking  control 
only  and  with  both  phase  and  tracking  controls  are  shown.  The 
ability  of  the  system  to  lock  in  a  nonoptimum,  off-boresight  con¬ 
dition  is  evident. 
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two-axis  focus  control  in  the  simulation,  in  effect  using  two  cylindrical  lens 
controls  that  focus  along  orthogonal  axes.  This  type  of  focusing  should  be 
particularly  useful  against  thermal  blooming  since  the  negative  thermal 
lens  is  much  stronger  perpendicular  to  the  wind  direction  than  parallel 
with  it. 

The  computer  simulation  results  have  shown  the  expected  10  msec 
convergence  time  for  the  tracking  controls  and  100  msec  for  the  focus. 

With  the  focus  control,  however,  we  observed  another  initially  unexpected 
result:  the  final  converged  level  can  depend  on  the  magnitude  of  the  initial 
defocus.  This  behavior  is  now  understood  to  be  a  generalized  manifestation 
of  the  ,!2mr  problem1'  encountered  with  deformable  mirror  COAT  systems. 
This  simulation  has  been  used  to  study  turbulence  and  blooming  compensa¬ 
tion;  results  of  these  calculations  are  discussed  in  Section  IV. 

C.  Thermal  Blooming  Generators 

Several  different  absorption  cells  have  been  built  during  the  course 
of  this  contract  and  each  one  has  been  used  in  a  variety  of  experimental 
arrangements.  Both  liquid  and  gaseous  media  have  been  used.  We  did  not 
observe  any  medium-dependent  differences  in  the  nature  of  blooming  or 
the  amount  of  COAT  correction  that  could  be  produced. 

1 .  Gas  Cells 

5 

Several  designs  were  evaluated  for  a  cell  to  produce 

7 

convection-dominated  thermal  blooming.  The  first  design  we  wmployed  is 
shown  in  Fig.  8.  The  design  uses  a  stationary  gas  cell  (Fig.  8(a))  and  a 
moving  mirror  arrangement  to  simulate  a  transverse  wind  (Fig.  8(b)). 
Experiments  in  which  the  cell  was  mounted  vertically  to  minimize  convec¬ 
tion  effects  gave  substantially  the  same  results  as  those  in  which  the  cell 

was  horizontal  (see  Section  III  of  this  report). 

5 

The  initial  cell  design  identified  SF^,  seeded  with  NO^  for  absorp¬ 
tion,  as  the  best  gas  for  producing  maximum  blooming  effects.  For  reasons 
that  are  still  not  apparent,  SF^  did  not  produce  the  expected  amount  of 
blooming.  We  theorize  that  the  large  number  of  V-T  states  in  SF^  pre¬ 
vents  rapid  transfer  of  the  energy  absorbed  by  the  NO^  into  thermal 
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Fig.  8.  Statis  gas  absorption  cell,  (a)  Photograph  of  cell. 

(b)  Schematic  of  "wind"  generation  apparatus.  The  moving 
mirror  causes  the  input  beam  to  move  relative  to  the  gas 
cell.  The  fixed  mirrors  bring  the  moving  beam  to  the 
back  of  the  moving  mirror  so  that  the  output  beam  is 
stationary. 


(translational)  energy  in  the  SF^.  The  next  choice  for  a  buffer  gas  was  Xe 
and  it  produced  blooming  according  to  our  predictions.  The  disadvantage 
of  Xe  is  its  much  higher  cost  compared  with  SF^. 

Although  the  static  gas  cell  in  Fig.  8  operates  satisfactorily,  it  has 
two  serious  drawbacks.  First,  it  cannot  be  used  with  slewing  beams  because 
of  limited  input  and  output  aperture  sizes.  More  important,  however,  the 
optical  arrangement  places  the  blooming  distortions  in  the  last  55%  of  the 
focused  propagation  path.  This  is  not  only  an  unrealistic  propagation 
scenario,  but  as  discussed  in  Section  III,  it  precludes  observing  any  sig¬ 
nificant  COAT  compensation  for  blooming.  To  correct  these  deficiencies, 

6  8 

we  built  a  flowing  gas  cell  9  as  shown  in  Fig.  9. 

The  flowing  gas  cell  can  produce  transverse  wind  velocities  up  to 
12  cm/sec  as  measured  by  a  hot-wire  anemometer.  The  fill  gases  are 
again  Xe  and  NO^.  No  turbulence  effects  are  observed  indicating  good 
laminar  flow.  This  cell  was  used  in  the  same  vertical  configuration  as 
the  static  gas  cell  to  minimize  buoyancy  effects  although  our  observations 
have  indicated  that  this  is  an  unnecessary  complication  of  the  optics.  The 
optical  arrangement  is  similar  to  Fig.  8(b)  except  that  the  focusing  lens 
can  be  placed  much  closer  to  the  cell  input  window  and  the  mirror  labeled 
1  or  2  in  Fig.  8(b)  is  stationary. 

We  experienced  some  difficulties  in  using  the  flowing  gas  cell 
mounted  as  shown  in  Fig.  9(b).  First,  the  blower  was  not  vacuum  tight 
so  that  small  quantities  of  NO^  slowly  leak  out  as  it  runs.  The  leak  is 
too  small  to  cause  a  safety  hazard,  but  it  is  annoying  because  of  the 
unpleasant  smell  of  NO^.  Second,  the  optics  were  cumbersome  and 
numerous  and  had  just  enough  residual  distortion  to  prevent  accurate 
measurement  of  the  beam  distortion  produced  by  blooming  alone  with  no 
COAT  correction.  Since  natural  convection  effects  are  not  noticeable 
for  our  experimental  conditions,  a  horizontal  mounting  arrangement  could 
be  used.  This  change  will  probably  eliminate  the  optics  problem,  but  will 
not  help  the  NO^  leak  problem.  Time  on  this  contract  did  not  permit  us  to 
try  out  a  horizontal  mounting,  but  we  expect  to  use  such  an  arrangement 
in  future  planned  experiments. 
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(a)  Schematic  drawing.  Sintered 
nickel  plates  are  used  for  dif¬ 
fusers  to  produce  laminar 
flow. 


Ml  1189 


(b)  Ceil  is  mounted  with  a  vertical  path  to 
minimize  buoyancy  effects. 


Fig.  9.  Flowing  gas  cell  for  blooming  studies. 
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2. 


Liquid  Cells 


All  of  the  thermal  blooming  studies  on  this  program  have 
been  scaled  to  simulate  conditions  encountered  by  focused  beams  propagating 
over  relatively  short,  low  altitude  paths  in  the  atmosphere.  The  scaling 
considerations,  discussed  in  detail  in  Section  III-A,  assume  that  the  heat 
transfer  is  dominated  by  forced  convection  (wind),  which  is  the  normal  con¬ 
dition  in  the  real  atmosphere.  As  long  as  this  assumption  is  satisfied,  and 

9 

if  we  are  not  concerned  with  effects  like  convective  cooling  that  occur  only 
at  high  altitudes,  there  are  no  fundamental  differences  between  liquid  and 
gaseous  blooming  media.  The  differences  that  have  been  observed^  occur 
only  for  collimated-beam  propagation  and  not  for  focused  beams  such  as  we 
are  considering. 

Liquid  media  have  two  distinct  advantages  over  gases  for  our  visible- 
wavelength  experiments.  First,  a  flowing-liquid  cell  is  simple  to  construct 
and  the  liquids  are  nontoxic  and  easy  to  handle.  Second,  much  stronger 
blooming  effects  can  be  produced  for  a  given  optical  power  and  absorption 
because  liquids  are  available  with  large  values  of  dn/dT,  the  temperature 
coefficient  of  the  refractive  index. 

As  part  of  the  Hughes  Research  Laboratories  IR&D  program,  we 
have  studied  COAT  compensation  for  thermal  blooming  in  liquids.  The 
results  obtained  are  included  with  this  contract  work  because  of  their 
applicability  to  the  questions  being  asked  on  this  program  and  because  the 
DARPA/RADC  COAT  system  was  used  in  the  measurements.  Two  liquid 
cells  have  been  used:  a  1  mm  thick  static  cell,  and  a  20  cm  long  flowing 
liquid  cell.  The  flowing  cell  is  shown  in  Fig.  10.  We  have  used  carbon 
tetrachloride  (CCl^)  methanol  for  liquids  with  iodine  added  to  absorb  the 
0.488  pm  laser  light.  Methanol  is  more  convenient  to  use  than  CCl^  in  the 
flowing  cell  because  the  very  high  value  of  dn/dT  in  CCl^  makes  it  hard  to 
produce  a  uniform,  turbulence-free  flow. 

D.  Turbulence  Phase  Screen 

Generating  artificial  turbulence  in  the  laboratory  is  not  difficult,  but 
doing  it  in  a  calibrated  fashion  is.  We  chose  to  use  a  single  phase  screen 
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4640-1 


Fig.  10.  Flowing  liquid  cell  for  thermal  blooming  studies. 

(a)  Schematic  diagram,  (b)  Photograph  of  20  cm  cell. 
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which  is  placed  just  ahead  of  the  lens  which  focuses  through  the  blooming 
cell.  The  phase  screen  is  fabricated  as  shown  in  Fig.  11(a)  by  sputtering 
patches  of  SiC>2  onto  a  fused  quartz  plate  through  a  mask.  The  mask  is  a 
thin  metal  plate  with  equal  size  holes  drilled  in  it  in  a  regular  pattern. 
Three  to  six  evaporations  of  about  1200  A  of  SiC>2  (\/4  at  4880  A)  are  made 
with  the  mask  in  a  different  position  each  time.  In  this  way  the  plate  is 
nearly  covered  with  a  random  patch  pattern  and  the  phase  variations  range 
from  zero  to  as  much  as  3X./2.  An  electronic  micrometer  scan  through  one 
of  these  phase  plates  is  shown  in  Fig.  11(b).  The  peak-to-peak  phase 
variation  is  about  3/4  wavelength  and  the  minimum  size  of  a  given  "turbule" 
is  about  0.  5  mm.  The  size  of  an  element  in  the  COAT  array  is  also  about 
0*  5  mm  at  the  location  where  the  plate  is  used.  To  produce  time-varying 
distortions,  the  plate  is  rotated  through  the  beam  as  shown  in  Fig.  11(b). 

This  type  of  turbulence  phase  screen  has  two  disadvantages.  First, 
the  turbulence  occurs  only  in  one  plane,  unlike  atmospheric  turbulence 
which  is  a  distributed  path  effect.  Second,  the  magnitude  of  the  distortions 
and  their  spatial  frequency  spectrum  are  quite  different  from  atmospheric 
turbulence.  In  particular,  we  have  found  that  our  phase  screens  contain 
very  little  tilt  or  beam  steering  distortion;  tilt  is  one  of  the  larger  compo- 

1  j 

nents  of  atmospheric  turbulence. 

The  artificial  nature  of  the  turbulence  is  not  important  for  the  pur¬ 
poses  of  this  program,  however.  We  have  already  demonstrated  the  ability 
of  COAT  to  correct  for  turbulence5'  12  and  we  have  shown  analytically^’  13 
that  turbulence  and  blooming  effects  can  be  considered  separately  and 
then  combined  later  to  give  an  effective  overall  reduction  in  beam  Strehl 
ratio.1  In  this  contract,  we  are  investigating  the  correction  for  turbulence¬ 
like  effects  in  the  presence  of  thermal  blooming. 

E.  Laboratory  Optical  Arrangement 

The  details  of  the  optical  layout  differ  for  each  experiment.  Each 
arrangement  has  several  common  features,  however,  and  in  functional 


*  Strehl  ratio  is  the  ratio  of  the  peak  beam  irradiance  to  the  diffraction- 
limited  irradiance. 
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.  11.  Phase  screen  for  generating  artificial  turbulence 

effects,  (a)  Evaporated  random  patches,  SiO^  on  a 
fused  quartz  substrate,  (b)  Calibration  scan  through 
an  artificial  turbulence  plate  using  an  electronic 
micrometer.  A  change  in  thickness  of  4880  R  will 
produce  a  one  wavelength  phase  shift  for  our  experi¬ 
ments .  (c)  Time  variations  produced  by  rotating  plate 
through  the  laser  beam. 


block  form  they  all  look  like  the  schematic  in  Fig.  12.  The  first  reducing 

telescope  recollimates  the  transmitting  array  to  a  size  usable  by  the  BAT 

actuators.  There  are  two  optical  paths  that  have  identical  optics  and  differ 

only  in  that  one  path  has  the  turbulence  and  blooming  distortions  in  it.  The 

second  undistorted  path  is  required  so  that  the  COAT  system  can  be  used  to 

1  2 

remove  random  step  phase  differences  introduced  by  the  nphasor  matrix.”  9 

When  the  photomultiplier  (PMT)  at  location  No.  2  is  connected  to  the 
COAT  servo  system,  a  diffraction-limited  beam  is  put  into  the  path  contain¬ 
ing  the  BAT  actuators  and  the  distortions.  The  Strehl  ratio  with  no  COAT 
correction  is  measured  this  way.  The  degree  of  COAT  correction  is  mea¬ 
sured  by  switching  the  PMT  at  location  No.  1  into  the  servo  system.  This 
comparison  procedure  is  possible  only  if,  when  the  blooming  and  turbulence 
distortions  are  removed,  there  is  no  change  if  the  beam  at  the  primary 
target  plane  as  the  PMT  choice  is  switched  from  No.  2  to  No.  1  and  back. 

As  mentioned  previously  and  as  discussed  further  in  Section  III,  this 
requirement  has  limited  our  ability  to  measure  the  effectiveness  of  focus 
control  alone  in  compensating  turbulence  and  blooming. 

For  all  the  experiments  discussed  in  Section  III,  we  will  use  the 
following  terminology:  (1)  nNO  COAT  correction"  means  PMT  No.  2  in 
Fig.  12  is  used;  (2)  "COAT  Corrected"  means  PMT  No.  1  is  used.  Occa¬ 
sionally,  we  may  refer  to  "local  loop  back,"  meaning  use  of  PMT  No.  2, 
or  "target  loop  lock,  "  meaning  use  of  PMT  No.  1.  In  all  cases,  the  beam 
profile  and  the  Strehl  ratio  are  measured  at  the  "primary  target  plane" 
indicated  in  Fig.  12. 
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Fig.  12.  Block  diagram  of  experimental  arrangement  used  in 
all  COAT  experiments. 
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III. 


EXPERIMENTAL  TECHNICAL  ACCOMPLISHMENTS 


The  experimental  work  on  this  program  consists  of  five  major  tasks: 
(1)  construction  of  a  gas  cell  thermal  blooming  simulator,  (2)  construction 
of  multidither  tracking  and  focus  controls,  (3)  measurement  of  COAT  com¬ 
pensation  for  thermal  blooming  with  single  and  multiple  glint  targets, 

(4)  measurement  of  COAT  compensation  for  turbulence  and  thermal  bloom¬ 
ing  simultaneously  with  complex  targets,  (5)  preliminary  study  of  target- 
signature/ speckle  effects  on  a  multidither  COAT  system.  The  last  three 

tasks  also  have  couterparts  in  the  analytical  studies  discussed  in  Section  V. 

6  8 

The  previous  section,  as  well  as  earlier  reports,  -  summarized  the  per¬ 
formance  of  the  hardware  that  was  built.  This  section  summarizes  the 
results  of  the  experimental  measurement  tasks. 

A  .  Thermal  Blooming  Compensation 

1 .  Scaling  Considerations 

The  thermal  blooming  experiments  on  this  program  are 
designed  to  increase  our  understanding  of  thermal  blooming  and  of  the 
limitations  of  multidither  COAT  systems  in  trying  to  remove  these  beam 
distortions.  The  experiments  are  thus  scaled  so  that  the  results  can  be 
extended  to  other  wavelengths,  ranges,  etc.  Since  blooming  effects  are 
most  pronounced  at  10.6  pm,  we  chose  to  scale  our  visible  wavelength 
laboratory  experiments  to  a  typical  10.6  pm  scenario.  In  order  to  scale  a 
10.  6  pm,  long-path  experiment  to  a  visible  wavelength  laboratory  experi¬ 
ment,  three  things  must  be  considered:  (1)  wind  speed  and  beam  size  to 
ensure  convection-dominated  heat  transfer  (the  only  case  of  interest  with 
actual  atmospheric  conditions),  (2)  total  absorption,  and  (3)  blooming 
strength  (transmitter  power). 

The  wind  speed  and  beam  size  are  related  to  the  medium  properties 
by  the  Peclet  number  defined  by 


2Vro 

Pe  =  (K/pC  ) 

r 


(1) 
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where  V  is  the  transverse  wind  velocity,  rQ  is  the  beam  radius  and  the 
medium  properties,  K,  p,  and  C  are  the  thermal  conductivity,  density, 
and  heat  capacity  at  constant  pressure.  The  quantity  D  =  (K/pC^)  is  also 
known  as  the  "diffusivity.  "  For  almost  all  atmospheric  conditions  of  inter¬ 
est,  heat  transfer  is  dominated  by  convection  rather  than  conduction.  This 
fact  is  equivalent  to  the  statement  that  negligible  energy  is  lost  by  a  mole¬ 
cule  in  the  time  required  for  the  molecule  to  move  across  the  optical  beam. 
Mathematically,  this  condition  can  be  expressed  as  Pg  »  1,  with  P  given 
by  eq.  (1).  Thus  for  a  given  medium,  there  is  a  minimum  value  for  the 
product  VrQ. 

To  ensure  that  the  blooming  distortions  as  a  function  of  propagation 
distance  are  the  same  in  the  laboratory  at  0.488  pm,  the  total  absorption, 
o-L,  and  the  Fresnel  number,  must  be  held  constant.  The  absorption 
coefficient  is  a  and  the  total  propagation  distance  in  the  blooming  medium  is 
L.  The  Fresnel  number  is  conveniently  defined  by 


kr 


(2) 


where  k  =  2tt/\  and  X  is  the  wavelength. 

5 

The  blooming  strength  is  measured  by  the  product  •  PT>  where 
P,j,  is  the  laser  power  at  the  entrance  to  the  blooming  medium  and  is 
defined  for  a  gaseous  medium  by 


Yo  (y  -  1 )  key  Le 

Trvr  V 
o  p 


-aL 


(3) 


The  quantities  y  , 
specific  heat  (y  = 
coefficient  for  the 


y,  and  p  are  the  molecular  polarizability  of  the  gas,  the 
C  /Cv),  and  the  static  gas  pressure.  The  absorption 
gas  is  a.  For  a  liquid,  a ^  is 
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where  the  liquid  density  is  p,  the  nominal  index  of  refraction  is  nQ,  and 
dn/dT  is  the  change  of  index  with  temperature. 

Wavelength  and  power  scaling  are  achieved  by  maintaining  a L  •  PT 
constant.  Table  3  presents  a  scaling  example  comparing  10.6  pm  propaga¬ 
tion  in  the  atmosphere  to  0.  488  pm  propagation  in  a  50  cm  cell  filled  with 
Xe  and  NO^.  -As  can  be  seen  from  Table  3,  scaling  to  a  450  kW,  10.  6  pm 
blooming  experiment  can  be  achieved  within  the  laboratory  using  reasonable 
experimental  parameters. 

As  mentioned  in  Section  II-C,  we  chose  xenon  seeded  with  NO.,  for 
our  experiments  with  gaseous  blooming  media.  A  partial  pressure  of  only 
3  to  4  Torr  of  NO?  is  required  to  achieve  50  to  70%  absorption  at  0.  488  pm. 

La 

The  cell  was  always  operated  with  a  total  pressure  of  one  atmosphere.  We 
initially  felt  that  meaningful,  scalable  experiments  could  be  accomplished 
only  with  a  gaseous  medium,  but  there  is  no  good  evidence  to  support  this 


Table  3.  Example  of  the  Scaling  Between  a  10.6  pm  Atmospheric 
Blooming  Experiment  and  a  0.488  pm  Laboratory  Gas 
Cell  Experiment 


Experiment 

Parameters 

10,6  pm  Atmospheric 
Experiment 

0.488  pm 

SF^+  NO^  Experiment 

Gas  Pressure 

1  atm 

1  atm 

0.  5 

0.  5 

L 

Z  km 

50  cm 

ro 

0.35m 

0.13  cm 

V 

1 0  m/  s 

6  cm/  s 

Pg  (eq.  23) 

4  x  105 

20 

°L 

2.66  x  10‘5 

1.14  x  102 

PT 

4.3  x 105  W 

0.1  W 

aL  PT 

11.4 

11.4 

T 1224 
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restriction  for  focused  beam  propagation.  As  data  in  the  next  section  show, 
similar  COAT  compensation  results  are  obtained  with  liquid  and  gaseous 
media  when  the  experimental  conditions  are  properly  scaled.  Since  strong 
blooming  is  much  easier  to  achieve  in  liquids,  there  is  little  motivation 
for  using  a  gas  blooming  cell  for  future  experiments. 

2.  Experimental  Results 

Several  experimental  measurements  of  COAT  compensation 
for  blooming  have  been  made.  The  principal  differences  between  each  pair 
of  experiments  are  the  ratio  of  the  length  of  the  blooming  medium  to  the 
focused  propagation  path  and  the  location  of  the  blooming  medium  within  the 
path. 

The  first  experiment  was  actually  part  of  the  Hughes  IR&D  program. 
The  experiment,  shown  schematically  in  Fig.  13,  used  a  thin  liquid  cell 
filled  with  CCl^  doped  with  iodine.  There  was  no  transverse  wind  in  this 
experiment.  The  basic  two-beam  arrangement  described  in  Section  II-E 
can  be  seen  in  Fig.  13.  When  the  liquid  cell  is  placed  as  close  as  possible 
to  lens  L^,  the  blooming  medium  is  in  the  collimated-beam,  near-field 
region  of  the  transmitting  array.  As  we  shall  see,  COAT  correction  has 
its  greatest  effect  on  near-field  blooming  distortions. 

With  the  cell  very  close  to  lens  L1  ,  the  results  shown  in  Fig.  14  are 
14  1 

obtained.  Both  the  boresight  irradiance  and  the  peak  irradiance  in  the 

beam  are  shown  as  a  function  of  power  into  the  cell.  The  theoretical  curve 

7  14 

shown  in  the  figure  is  obtained  from  a  simplified  theory  9  that  compares 
the  uncorrected  thermally-induced  beam  spreading  with  the  diffraction 
angle  of  a  collimated  beam.  The  value  of  the  parameter  A  is  obtained  by 
fitting  the  theory  to  the  experiment  at  one  power  level.  In  this  experiment, 
the  transmitted  beam  is  not  focused  to  produce  the  far-field  irradiance  pat¬ 
tern,  but  instead  expands  by  diffraction  to  produce  the  far-field  pattern. 

The  data  in  Fig.  14  indicate  that  the  COAT  system  increases  the 
optimum  transmitter  power*  by  a  factor  of  2  to  3  and  increases  the  peak 
target  irradiance  by  up  to  a  factor  of  4.  More  detailed  data  on  beam  pro- 

7 

files  have  been  presented  elsewhere.  We  hypothesize  at  this  point  that  the 


*The  transmitted  power  that  gives  maximum  target  irradiance. 
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ELECTRONICS 

Fig.  13.  Optical  paths  for  static  liquid  cell  blooming  experiment. 

BS:  beam  splitter.  M:  mirror.  PMT:  photomultiplier. 
Lens  LQ  together  with  either  lens  Lj  or  L2  forms  a  lOOx 
reducing  telescope . 
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Fig.  14.  Thermal  blooming  compensation  with  a  thin,  static-liquid 
blooming  cell  located  in  the  near-field  of  the  transmitter. 
The  dashed  curve  is  theoretical  (see  references  7  and  14) 
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COAT  system  is  this  effective  only  because  the  blooming  distortions  occur 

in  the  near-field  of  the  transmitter  and  that  as  the  blooming  occurs  farther 

from  the  transmitting  aperture,  the  COAT  phase  correction  will  be  less 

7 

effective.  (In  fact,  a  simple  analysis  can  be  used  to  show  that  correction 
in  this  experiment  is  being  limited  mainly  by  the  finite  number  of  corrector 
elements  and  not  by  a  fundamental  inability  of  COAT  to  remove  the 
distortions.  ) 

To  test  this  hypothesis,  we  repeated  the  measurements  of  Fig.  14 
with  the  cell  at  different  distances  from  lens  L.  .  The  results  are  shown  in 

z  1 

Fig.  15.  The  quantity  ZR  =  DT/4\  is  the  classical  Rayleigh  distance, 
which  is  the  distance  at  which  the  transition  from  near-field  to  far-field 
occurs.  For  our  experiment,  ZR  =  9  mm.  As  expected,  the  degree  of 
COAT  correction  is  reduced  significantly  as  the  blooming  cell  is  moved 
toward  the  far-field  of  the  transmitting  aperture  (d/zR  £  !)•  In  fact,  it 
appears  that  significant  correction  occurs  only  for  distortions  located  well 
within  the  first  30%  of  the  propagation  distance  from  near-field  to  far-field 
since  no  appreciable  correction  is  observed  for  d/ Z^  >  0.28. 

The  next  experiment  employed  the  static  gas  cell  in  the  arrangement 
shown  in  Fig.  16.  In  this  case,  the  laser  beam  is  actually  focused  through 
the  cell.  The  data  obtained  with  this  arrangement  are  shown  in  Figs.  17 
and  18.  The  photographs  in  Fig.  17  are  black-and-white  photographs  from 
a  level-quantized  color  TV  display^  with  1.  5  dB  per  color  level. 

There  are  five  principal  features  of  these  data.  First,  the  expected 
crescent-shaped  beam  distortion  is  observed  in  Fig.  17(c)  and  the  beam  is 
shifted  into  the  wind.  Second,  when  the  COAT  correction  is  applied,  the 
beam  maximum  is  maintained  on  the  horesight  axis  where  the  glint  is 
located.  Third,  there  is  no  improvement  in  peak  irradiance  when  the  COAT 
correction  is  applied,  except  for  a  small  increase  in  the  power  through  an 
aperture  whose  width  equals  the  full-width  at  half-maximum  of  the  unbloomed 
beam.  Fourth,  the  system  convergence  time  was  unaffected  by  blooming 
level,  being  1.  5  to  2.0  msec  at  all  power  levels.  Finally,  the  convergence 
is  stable  at  all  power  levels. 
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Fig.  16.  Schematic  of  laboratory  optics.  Path  No.  1  through  the 

blooming  cell  and  path  No.  2  outside  the  cell  are  optically 
identical  when  the  cell  is  not  filled. 
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Fig.  17.  Forced  convection  thermal  blooming.  Each  shading 
level  differs  from  the  adjacent  one  by  1 . 5  dB . 
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8  15 

We  had  speculated  earlier  '  that  the  absence  of  correction  in  this 
experiment  may  have  been  caused  by  buoyancy  effects  or  by  too  low  wind 
velocities.  We  now  believe  that  we  see  no  correction  because  the  blooming 
medium  is  in  the  last  55%  of  the  propagation  path.  Strong  phase  distortions 
in  this  region  require  amplitude  changes,  not  phase  corrections,  at  the 
transmitting  aperture.  The  servo  response  speed  may  also  play  some  role 
(see  discussion  below),  but  this  is  not  clear  now.  However,  the  lack  of 
COAT  correction  in  data  in  Fig.  18  is  entirely  consistent  with  the  data  in 
Fig.  15. 

One  experimental  run  has  been  made  with  the  flowing  gas  cell.  The 
BAT  tracking  controls  were  included  in  these  tests.  The  experimental 
arrangement  was  similar  to  those  shown  in  Fig.  16,  except  that  lens  f^ 
(focusing  through  the  cell)  was  moved  and  its  focal  length  selected  so  that 
the  blooming  occurs  in  the  last  81%  of  the  propagation  path.  The  observed 
improvement  in  peak  irradiance  is  a  factor  of  1. 35  or  35%.  This  improve¬ 
ment  was  produced  by  the  phase  controls  alone;  the  tracking  controls  either 
alone  or  when  added  to  the  phase  controls,  produced  no  observable  change 
in  the  peak  target  irradiance. 

As  mentioned  previously,  we  were  not  able  to  determine  the  com¬ 
pensation  of  thermal  blooming  using  the  dithered  focus  control  alone 
because  of  spherical  distortion  in  the  focus  actuator.  We  did  try  it  as  an 
addition  to  the  18-channels  of  phase  control,  however,  but  we  observed  no 

change  in  the  peak  target  irradiance  when  the  focus  control  loop  was  closed. 

16-18 

This  result  is  consistent  with  other  work  “  that  has  demonstrated  little 
or  no  improvement  with  spherical  focus  control  alone,  but  significant  cor¬ 
rection  with  coma  and  astigmatism  control. 

This  improvement  is  more  encouraging,  but  is  still  less  than  the 

16  18  20 

factors  of  1.75  to  2.0  reported  by  others  *  ’  The  comparison  is  unfair, 

however,  since  the  work  in  the  referenced  literature  used  slewing  beams 
(which  effectively  eliminates  blooming  near  the  focal  plane).  Our  experi¬ 
ments  had  no  slewing,  and  the  distortion  was  in  the  last  81%  of  the  path, 
thus  presenting  a  much  harder  task  for  the  COAT  system.  In  fact,  referring 
again  to  the  data  in  Fig.  15,  we  observed  a  factor  of  1 . 7  improvement  in 
peak  target  irradiance  when  the  blooming  occurred  roughly  in  the  last  83% 
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of  the  propagation  path  (compare  the  uncorrected  curve  to  the  corrected 

curve,  d/Z„  =  0.17,  in  Fig.  15).  Actually,  in  the  experiment  leading  to 

Fig.  15,  the  blooming  occurred  only  from  d/Z^  =  0.17tod/Z^  =  0.28 

rather  than  over  the  last  83%  of  the  "focused"  path  (d/ ZR  =  0.  17  to 

d/Z  =  1.0).  Again,  this  is  a  less  severe  case  as  far  as  COAT  correction 
R. 

is  concerned  than  the  experiment  giving  Fig.  19. 

The  final  blooming  experiment  used  a  20  cm  long  flowing  liquid  cell. 
Again,  an  optics  limitation  prevented  us  from  placing  the  entire  focused 
beam  path  within  the  blooming  medium.  The  experimental  arrangement  is 
similar  to  that  shown  in  Fig.  12.  Two  experiments  were  performed.  In 
the  first  experiment,  the  blooming  medium  was  in  the  last  72%  of  the 
propagation  path.  The  results  of  these  experiments  are  shown  in  Fig.  20. 
Not  unexpectedly,  the  data  are  very  similar  to  those  in  Fig.  19,  and  we 
observe  a  factor  of  1.33  improvement  in  peak  target  irradiance.  This 
result  gives  us  confidence  that  there  are  no  significant  differences  between 
focused-beam  COAT  compensation  experiments  performed  with  liquid  and 
gaseous  media. 

In  the  second  experiment,  the  liquid  cell  was  moved  to  the  first  72% 
of  the  path,  producing  the  data  in  Fig.  21.  We  now  see  a  factor  of  1 .  5 
increase  in  the  optimum  peak  target  irradiance.  A  factor  of  2  increase 
from  the  uncorrected  to  the  corrected  case  is  seen  at  the  highest  power 


level.  This  arrangement  produces  blooming  and  compensation  results 

along  the  propagation  path  that  are  similar  to  that  occurring  in  a  slewed 

^  ,.  16,19,20  _ 

beam  when  the  entire  propagation  path  is  in  the  medium.  Experi¬ 

ments  using  stationary  and  slewed  beams  that  are  focused  through  the 
entire  blooming  medium  will  be  performed  as  part  of  follow-on  work  to 
this  contract  (contracts  F30602-76-C-0021  and  F30602-76-C-0022). 


3.  Conclusions 

We  now  feel  that  our  closed-loop,  real-time  phase  compensa¬ 
tion  results  are  completely  consistent  with  other  results  that  were  obtained 

16  19  20 

using  predictive,  open-loop  methods.  *  1  A  factor  of  about  1.  5  to 

2.0  improvement  in  peak  target  irradiance  at  optimum  transmitter  power  is 
possible  with  multidither  COAT  dependong  on  the  slewing  rate.  Correction 
factors  of  2  or  more  may  be  possible  at  higher  transmitter  powers,  but  a 
lower  peak  target  irradiance  is  produced. 
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Fig.  19.  Thermal  blooming  compensation  data  obtained  with  flowing  gas 
cell.  Also  shown  is  the  focused  propagation  path  relative  to 
the  blooming  medium. 
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Fig.  20.  COAT  thermal  blooming  compensation  with  a  flowing  liquid 

medium.  The  medium  is  in  the  last  72%  of  the  path.  Results 
are  similar  to  those  in  Fig.  19. 
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Fig.  21.  COAT  correction  with  a  flowing  liquid  medium 
that  is  in  the  first  72%  of  the  propagation  path. 
Significantly  better  performance  is  observed 
than  in  the  experiment  described  in  Fig.  20. 
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Care  must  be  exercised  when  comparing  compensation  results 

from  different  experiments,  since  our  data  indicate  that  only  the  blooming 

distortions  which  occur  in  the  first  20  to  30%  of  the  focused  propagation  path 

can  be  removed  by  correcting  the  transmitter  phase.  This  suggests  that 

any  method  that  can  minimize  the  blooming  distortions  close  to  the  focal 

plane  will  be  effective  not  only  in  increasing  the  target  irradiance,  but  also 

in  making  the  remaining  distortions  more  amenable  to  COAT  correction. 

Three  techniques  come  immediately  to  mind:  (1)  transmitter  intensity 

tailoring  or  control,  “  (2)  use  of  slewed  beams,  (3)  use  of  pulsed 

25 

lasers  with  repetition  rates  low  enough  to  allow  only  2  to  3  pulses  during 
a  wind-driven  medium  transit  time  across  the  beam.  When  comparing 
improvement  factors  in  different  experiments,  we  must  also  be  careful  to 
specify  whether  we  are  comparing  the  maximum  achievable  target  irradiance 
with  and  without  correction  or  the  improvement  at  a  fixed  transmitter  power 
level.  The  use  of  slewing  beams  also  greatly  affects  the  comparisons. 

In  comparing  our  blooming  experiments  to  COAT  system  tests  done 

3  4  5 

with  no  distortions  9  or  with  atmospheric  turbulence  in  the  path,  we  have 

found  no  differences  in  system  convergence  time,  multiple  glint  discrimina¬ 
tion,  or  tracking  ability.  We  checked  these  quantities  carefully  in  the  last 
two  experiments  described  above.  We  also  noted  that  the  same  blooming 
compensation  was  observed  with  a  diffuse- surface  sphere  for  a  target 
instead  of  a  small  point-glint  target. 

Although  the  thermal  blooming  distortions  do  not  affect  the  COAT 
servo  response  time,  we  are  concerned  about  what  the  response  time  should 
be  for  optimum  blooming  compensation.  Analysis  of  this  problem  is  greatly 
complicated  by  the  nonlinear  distortion  medium.  At  this  point,  let  us  just 
state  the  problem;  we  will  return  to  it  again  in  Section  IV. 

Since  the  COAT  system-plus -blooming  medium  is  a  multiple-loop, 
nonlinear  system,  it  is  not  unlikely  that  multiple  solutions  exist  for  the 
transmitter  phase  front  that  satisfies  the  servo  error  criterion:  near-zero 
dither  fundamental  frequency  in  the  signal  reflected  from  the  target.  The 
problem  is  to  choose  the  "path"  in  "correction  space"  that  ensures  the  sys¬ 
tem  stops  on  the  true  maximum  target  irradiance,  not  on  a  "local"  maximum. 
It  is  important  to  point  out  at  this  juncture  that  we  are  only  hypothesizing 
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the  existence  of  subsidiary  maxima  in  the  COAT  convergence  process  with 
a  nonlinear  distortion;  there  is  no  experimental  or  analytical  evidence  yet 
to  support  this  hypothesis. 

If  the  COAT  system  convergence  time  is  faster  than  the  response 
time  of  the  nonlinear  medium,  it  will  not  sample  the  true  slope  of  the  target- 
irradiance  versus  phase-error  curve  properly.  That  is,  the  instantaneous 
distortion  will  be  completely  removed.  If  instead  of  this  fast  response, 
the  COAT  system  was  slower  than  the  medium,  it  would  act  like  a  "step, 
look,  decide,  step"  servo  system.  In  this  type  of  system,  a  change  in  the 
transmitter  phase  front  is  made,  the  medium  responds  to  the  change,  and  a 
decision  is  made  whether  the  change  increased  or  decreased  the  target 
irradiance.  If  the  irradiance  increases,  the  direction  of  the  next  change  is 
the  same;  if  not,  the  direction  reverses.  This  is  effectively  the  type  of 
system  studied  by  others.1^’  The  amplitude  of  the  test  changes  (the 

dither  amplitude  in  a  multidither  system)  will  also  be  important  in  this  type 
of  system  if  local  maxima  are  to  be  avoided. 

At  this  point,  we  speculate  that  it  may  be  important  to  allow  the 
medium  to  fully  respond  to  the  COAT-induced  phase  changes  and  possibly 
even  to  the  dither  frequencies  if  a  multidither  COAT  system  is  to  converge 
to  the  global  optimum  transmitter  phase.  This  requirement  means  dither 
frequencies  well  below  the  present  8  to  32  kHz  band  used  in  the  DARPA/ 
RADC  system.  The  problem  with  such  a  slow  system  is  that  it  is  totally 
inadequate  for  turbulence  compensation.  A  possible  solution  to  this 
dilemma,  if  our  speculation  is  correct,  would  be  the  use  of  some  sort  of 
predictive  correction  so  that  the  COAT  system  starts  close  enough  to  the 
correct  phasefront  that  local,  sub-optimal  maxima  are  avoided.  A  second 
solution  may  lie  in  the  existence  of  turbulence  in  the  optical  path.  The 
random  nature  of  turbulence,  and  a  fast  COAT  system's  action  to  correct 
for  it,  may  be  enough  to  ensure  that  a  broad  region  of  "correction  space" 
is  sampled  so  that  the  global  optimum  phase  correction  is  always  found. 

It  is  important  to  point  out  that,  all  of  the  above  discussion  aside, 
we  have  never  observed  any  servo  instability  in  our  experimental  bloom¬ 
ing  compensation  measurements.  In  each  experiment,  the  servo  gain 
(response  time)  is  optimized  with  no  blooming  or  turbulence  distortions 
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present.  The  distortions  are  then  added  and  the  laser  power  is  increased 
(for  blooming).  In  no  instance  did  the  COAT  servo  become  unstable. 

The  blooming  compensation  may  have  been  less  than  anticipated,  but  a  stable 
phasefront  correction  was  always  achieved.  This  stable  behavior  is  not 
unexpected  for  intensity-maximization  systems  such  as  multidither,  out¬ 
going  wave  COAT,  but  may  not  be  characteristic  of  return-wave  systems 

26  27 

such  as  multidither,  return-wave  or  phase-conjugate  or  "TRIM-COAT11* 
systems. 

B.  Turbulence  Plus  Blooming  Compensation 

One  of  the  goals  of  this  contract,  and  an  important  consideration 
for  any  operational  COAT  system,  is  the  determination  of  whether  thermal 
blooming  will  diminish  a  COAT  system1  s  ability  to  correct  for  turbulence. 

This  section  presents  the  results  of  our  studies  of  this  problem. 

1 .  Artificial  Turbulence 

The  artificial  turbulence  phase  screen  was  produced  as  dis¬ 
cussed  in  Section  II-D  and  located  as  shown  in  Fig.  12.  The  distortions 
proved  to  be  slightly  too  strong  for  this  COAT  system.  That  is,  more  and 
smaller  elements  are  required  to  fully  compensate  for  the  distortions. 

With  no  COAT  correction,  a  time-averaged  (10  sec)  Strehl  ratio  of  0.22  was 
observed.  The  COAT  system  was  able  to  correct  this  beam  back  to  a  Strehl 
ratio  of  0.60.  The  phase  plate  was  rotated  to  produce  time-varying  turbu¬ 
lence.  A  plot  of  the  control  voltage  for  one  of  the  transmitting  array  ele¬ 
ments  is  shown  in  Fig.  22.  Significant  distortions  (as  measured  by  the 
corrections  for  them)  occur  at  frequencies  up  to  30  Hz.  This  is  in  reasonably 

good  agreement  with  the  50  Hz  spectral  width  observed  with  real  atmospheric 

5 

turbulence  with  this  COAT  system. 

Using  a  somewhat  less  severe  phase  screen,  we  have  compared  the 
use  of  phase  control  to  the  use  of  tracking  control  only.  This  phase  screen 
not  only  had  lower  phase  distortions,  but  the  quartz  plate  was  mounted  in 
such  a  way  that  significantly  more  steering  error  was  introduced  than  with 


^TRIM-COAT  is  the  acronym  applied  to  broad-band  image-compensation 
(I-COAT)  systems  that  are  used  to  correct  the  phase  of  a  transmitted, 
coherent  beam. 
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Fig.  22.  Artificial  turbulence  frequency  spectra.  The 
observed  signal  is  the  control  signal  from 
COAT  channel  No.  5.  (a)  Total  spectrum  to 

150  Hz,  (b)  Expanded  scale  spectrum  from  dc 
to  40  Hz.  Note  the  different  vertical  scales  for 
the  two  curve. 
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the  phase  screen  shown  in  Fig.  12  and  used  in  the  turbulence  plus  blooming 
measurements  discussed  below. 

We  have  looked  at  four  cases:  (1)  no  COAT  control,  (2)  phase  con¬ 
trol  only,  (3)  tracking  control  only,  (4)  phase  and  tracking  control  together. 
The  results  of  the  measurements  are  shown  in  Table  4.  As  can  be  seen 
from  the  data  in  the  table,  the  phase  and  steering  distortions  in  this  artifi¬ 
cial  turbulence  screen  have  roughly  equal  contributions  to  the  uncompensated 
distortion. 

2.  Turbulence  Plus  Blooming 

The  turbulence  phase  screen  used  in  these  tests  was  the 
stronger  of  the  two  discussed  above.  The  thermal  blooming  scenario 
chosen  was  the  same  as  that  producing  the  data  in  Fig.  21:  a  flowing  liquid 
medium  in  the  first  72%  of  the  propagation  path.  The  observed  compensa¬ 
tion  performance  is  shown  in  Fig.  23.  The  peak  irradiance  for  turbulence 
plus  blooming  is  almost  as  good  as  that  achieved  in  the  absence  of  turbu¬ 
lence.  The  difference  is  caused  by  the  inability  of  the  COAT  system  to  com¬ 
pletely  remove  the  turbulence.  This  can  be  seen  by  noting  that  initial  slope 
of  the  £  0,  COAT-corr ected  curve  in  Fig.  22  is  lower  than  the  slope  of 
the  free-space  line  and  by  noting  that  the  ratio  of  the  irradiances  for  the 
two  corrected  cases  at  any  power  level  is  very  close  to  0.60,  the  no-blooming, 
corrected  Strehl  ratio. 


Table  4.  COAT  Compensation  for  Artificial  Turbulence 


Case 

Peak  Target 
Irradiance, 
arb.  units 

Strehl  Ratio,  * 

S 

NO  COAT 

0.  39 

0.65 

COAT  phase  control 

0.  52 

0.87 

COAT  tracking  control 

0.  52 

0.87 

| 

Phase  and  tracking  control 

0.  57 

0.95 

*S  =  1  is  defined  for  no  turbulence,  COAT  phase  control. 
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Fig.  23.  COAT  correction  for  turbulence JLn  the  presence 
of  blooming  distortions.  The  =  0  curves  are 
the  same  as  in  Fig.  21  . 
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3. 


Conclusions 


The  data  presented  above  indicate  that  nearly  100%  of 
turbulence  distortions  can  be  removed  by  a  COAT  System  (given  a  sufficient 
number  of  COAT  channels)  even  though  the  absolute  Strehl  ratio  is  severely 
reduced  by  blooming.  The  nonlinear  blooming  medium  does  not  limit  the 
COAT  system's  ability  to  correct  for  turbulence.  This  conclusion  is  also 
supported  by  the  computer  simulation  data  discussed  in  Section  IV. 

C.  Complex  Target  Effects 

Contract  Amendment  No.  2  redirected  some  of  the  contract  work  to 
provide  time  and  funds  for  some  initial  studies  of  the  interaction  of  a  multi¬ 
dither  COAT  servo  system  with  backscatter  from  an  extended,  semidiffuse 

target.  Some  of  the  results  of  these  studies  were  presented  at  the  1975 

28 

Conference  on  Laser  Engineering  and  Applications.  The  interest  in 

pursuing  these  studies  was  initially  motivated  by  the  work  of  Ogrodnik  and 

29 

Gurski.  A  discussion  of  the  results  of  these  studies  can  also  be  found 
in  Ref.  6. 

1 .  Problem  Statement 

The  physical  problem  known  as  ’’speckle  noise”  can  be  stated 
briefly.  When  a  coherent  laser  beam  illuminates  a  target,  some  of  the 
energy  is  scattered  back  in  the  direction  of  the  transmitted  beam.  With  all 
except  single  glint  targets,  the  backscattered  radiation  produces  a  random 
intensity  pattern  (a  ’’speckle  pattern”)  at  a  receiver,  which  is  located  close 
to  the  transmitter  for  cases  of  interest  to  us.  When  the  appropriate  condi¬ 
tions  exist,  any  movement  of  this  speckle  pattern  relative  to  the  receiver 
will  produce  an  amplitude  modulation  of  the  received  signal.  Since  a  multi¬ 
dither  COAT  system  receiver  senses  amplitude  modulations  on  the  back- 
scattered  beam  that  are  produced  by  dithering  the  transmitted  beam  phase, 
large  amplitude  spurious  amplitude  modulations  in  the  receiver  may  swamp 
the  desired  modulations  and  thus  interfere  with  the  system  operation. 

Doppler  shifts  produced  by  rapidly  rotating  targets  may  also  produce 
false  signals,  but  these  spurious  signals  are  expected  to  be  less  important 
than  amplitude  modulation  effects  for  the  target  scenarios  of  interest.  This 
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is  true  because  the  most  likely  glint  point  on  a  rapidly  rotating  target  is  one 
normal  to  the  beam  incidence  direction;  this  point  will  have  nearly  zero 
transverse  doppler  shift.  Longitudinal  doppler  shifts  are  very  narrow  band 
and  thus  will  have  minimal  effect  on  a  multiple- channel  COAT  system. 

Thus  for  COAT  applications,  it  is  important  to  determine  three 
things  about  target-signature  effects.  First,  what  kind  of  modulations  inter¬ 
fere  with  the  proper  operation  of  a  multidither  COAT  system?  The  interest¬ 
ing  parameters  of  the  modulation  are  amplitude,  frequency  spectrum,  and 
power  spectral  density.  Second,  what  kind  of  receiver  modulations  are  pro¬ 
duced  by  backscatter  from  real  targets  in  expected  operational  scenarios? 
That  is,  what  scenarios  will  produce  the  modulations  that  can  affect  the 
COAT  system?  Third,  can  a  multidither  COAT  system  be  designed  that  is 
insensitive  to  speckle-induced  receiver  modulation  and  if  so,  what  are  the 
critical  design  parameters?  A  follow-on  contract.  No.  F30602-76-C-0021 , 
to  the  current  contract  is  designed  to  produce  definitive  answers  to  all  these 
questions.  Some  preliminary  studies  of  the  magnitude  of  the  effects  were 
obtained  on  this  contract,  however. 

2 .  Experimental  Results 

For  a  first  experiment,  we  chose  the  simple  extended  moving 
target  shown  in  Fig.  24(a).  The  target  is  a  cylindrical  scotchlite  strip  that 
could  be  rotated  at  different  rates.  The  target  is  larger  than  the  COAT- 
formed  array  pattern  in  both  dimensions  (parallel  with  and  perpendicular  to 
the  rotation  axis),  but  exceeds  the  element  pattern  diameter  only  in  one 
dimension.  The  curvature  of  the  scotchlite  provides  a  localized  highlight'1' 
on  which  the  COAT  system  can  lock.  The  laboratory  experimental  arrange¬ 
ment  used  is  shown  in  Fig.  12.  No  turbulence  or  blooming  was  introduced 
in  the  paths  and  the  BAT  controls  were  not  used.  The  receiver/transmitter 
aperture  ratio  was  chosen  to  be  0.  66. 

♦It  is  erroneous  to  say  that  a  multidither  COAT  system  requires  a  ’’glint" 
for  convergence.  It  will  not  converge  on  a  flat  plate  (at  least  the  systems 
considered  here  will  not),  but  any  highlight  a  change  in  reflectivity  with 
beam  position  on  the  target)  will  suffice  to  provide  the  required  reference. 
We  have  established  this  fact  experimentally  with  several  target  shapes. 


65 


Fig.  24.  COAT  receiver  power  spectrum  with  extended,  moving 

targets,  (a)  Extended  target  model,  (b)  Received  signal 
spectrum  with  stationary  target,  (c)  Received  signal 
spectrum  with  moving  (rotating)  target. 
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The  experiment  consisted  of  observing  the  peak  target  irradiance  and 
the  COAT  receiver  signal  frequency  spectrum  when  the  target  was  stationary 
and  moving.  Figures  24(b)  and  24(c)  shows  the  observed  COAT  photomulti¬ 
plier  receiver  power  spectra  for  a  typical  case.  The  spectrum  analyzer 
bandwidth  was  200  Hz.  Each  photograph  consists  of  five  2-min  scans  over 
a  0  to  50  kHz  range.  Figure  24(b)  is  the  observed  spectrum  when  the  target 
is  stationary  and  Fig.  24(c)  is  the  spectrum  when  the  target  rotates.  The 
target  rotation  speed  was  chosen  to  give  the  largest  observed  change  in  the 
spectrum  over  the  dither  band.  For  the  stationary  target,  the  noise  level 
is  low  and  the  dither  signals  are  clearly  evident  above  the  noise.  When  the 
target  is  moving,  the  noise  level  increases  significantly  and  most  of  the 
dither  signals  are  no  longer  evident.  The  peak  target  irradiance  was  the 
same,  however,  whether  the  target  was  stationary  or  moving.  For  this 
simple  target,  if  the  COAT  system  could  converge  on  the  target  when  it  was 
stationary,  the  system  could  converge  equally  well  when  the  target  was 
moving.  For  every  target  and  target  motion  rate  that  we  have  tried  to  date, 
this  same  observation  has  been  made:  no  noticeable  effects  when  the  target 
moves. 

The  observation  of  no  effect  on  the  COAT  system  performance  is 
contrary  to  a  statement  made  in  the  third  management  report.  The  very 
earliest  experiments  did  seem  to  show  some  effect.  What  we  found,  how¬ 
ever,  was  that  very  large  extraneous  dither  signals  were  being  introduced 
into  the  COAT  receiver  by  ground  loops,  not  by  target  speckle  effects. 

Since  the  ground  loops  have  been  eliminated,  we  have  not  been  able  to  find 
a  target  which  ’’confused11  the  COAT  system  solely  because  of  its  motion. 

Since  we  have  seen  no  deleterious  COAT-system/ target-modulation 
effects,  a  natural  question  is.  What  kind  of  spurious  signals  can  reduce  the 
system  performance?  We  thus  undertook  an  experiment  whose  goal  was  to 
generate  large  amplitude  spurious  amplitude -modulated  signals  capable  of 
overloading  the  dither  servo  channels.  In  this  experiment,  the  target  was  a 
single  glint  and  a  large  receiver  aperture  was  used.  Spurious  amplitude 
modulations  were  artificially  introduced  into  the  COAT  system  by  putting  an 
acousto-optic  modulator  (AOM)  in  the  optical  beam  path,  as  shown  in  Fig.  25. 
The  AOM  amplitude  modulates  the  transmitted  beam  at  a  frequency  set  by 
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Fig.  25. 


Experimental  setup  to  determine  effect  of  spurious  amplitude 
modulations  in  a  COAT  receiver. 
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the  voltage- controlled -oscillator  (VCO).  Modulation  depths  of  64%  were 
observed  over  the  dither  frequency  band  (8  to  32  kHz)  with  a  triangular- 
wave  drive  to  the  VCO. 

Manually  tuning  the  VCO  through  the  dither  passband  sequentially 
overloaded  the  feedback  loop  at  each  dither  frequency  and  caused  a  10% 
reduction  in  peak  target  irradiance  (a  loss  of  proper  phasing  in  1  of  the  18 
COAT  channels).  Spectrum  analyzer  measurements  indicated  that  the 
spurious  signal  level  of  100  mV  generated  in  this  manner  was  much  larger 
than  the  3  mV  dither  amplitudes.  This  initial  result  establishes  the  fact 
that  external  modulation  applied  to  the  COAT  beam  could  overload  the  con¬ 
trol  circuitry  in  any  single  channel  and  cause  the  loss  of  that  channel  in 
forming  the  converged  beam. 

To  roughly  simulate  the  effect  of  worst-case  broadband  speckle 
modulations,  the  VCO  was  swept  over  a  selectable  bandwidth  with  sweep 
periods  of  1  msec  and  0.  5  msec  (shorter  time  than  the  COAT  system  con¬ 
vergence  time  of  1.5  to  2.0  msec).  During  this  series  of  tests,  the  COAT 
receiver  frequency  spectrum  was  recorded  and  the  power  level  in  the  con¬ 
verged  beam  was  measured  with  the  sweep  voltage  turned  on  and  off.  The 
center  frequency,  f  ,  was  selected  for  maximum  effect  on  the  COAT  system. 
A  graph  of  the  number  of  active  elements  lost  versus  modulation  bandwidth 
is  shown  in  Fig.  26.  The  sweep  bandwidths  in  Fig.  26  range  from  1  kHz  to 

30  kHz  (f  =  f  ±  BW/2),  but  the  ratio  of  the  peak  target  irradiance  with  the 
o 

spurious  modulation  to  that  without  it  (pmocj/Po^  is  nearly  constant.  A 

value  of  P  ,/P  =  0.  69  =  (15/18)^  corresponds  to  a  loss  of  proper  phas- 

mod  o  ^ 

ing  in  three  channels  of  the  average  and  (14/18)  =  0.60  corresponds  to  a 

loss  of  4  channels.  Similar  results  are  obtained  when  the  sweep  period  is 
reduced  to  0.  5  nsec. 

The  power  spectrum  of  the  triangular-wave  swept  modulation  used 
here  consists  of  a  sum  of  all  the  harmonics  of  the  sweep  frequency  spaced 
about  the  center  frequency: 


P  =  >  a  exp  [ i2Tr  (f  -  nf  )] 

m  /  j  n  v  o  sweep  J 

n=-oo 


(4) 
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PEAK  TARGET  IRRADIANCE,  NORMALIZED  UNITS 


MODULATION  BANDWIDTH,  kHz 

Fig.  26.  Effect  of  acousto-optic  modulation  on  DARPA/RADC  COAT 
system  as  a  function  of  total  modulation  bandwidth. 
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where  the  a  are  determined  in  part  by  the  sweep  bandwidth.  Although  this 
n 

type  of  power  spectrum  will  probably  not  be  characteristic  of  any  target- 
induced  speckle  spectrum,  it  does  have  several  properties  which  a  speckle- 
induced  receiver  spectrum  should  also  have.  First,  the  bandwidth  of  the 
extraneous  receiver  signals  is  limited,  but  will  almost  certainly  exceed 
35  kHz  in  cases  of  most  interest.  Second,  the  product  of  peak  power 
spectral  amplitude  and  bandwidth  is  constant.  What  this  means  is  that  the 
receiver  can  see  a  100%  average  modulation -depth  signal  at  only  a  single 
frequency;  as  the  modulation  bandwidth  increases,  the  average  receiver 
voltage  modulation  depth  (not  the  instantaneous  depth  or  amplitude)  goes 
down  as  the  reciprocal  of  the  square  root  of  the  bandwidth. 

Determination  of  the  actual  spectrum  of  a  real  speckle-induced 

receiver  signal  is  beyond  the  scope  of  this  contract.  As  a  worst  case, 

however,  consider  a  spectrum  of  equal  amplitude  sinusoidal  frequencies, 

spaced  by  100  Hz  over  a  35  kHz  bandwidth,*  which  is  close  to  a  worst-case 

situation.  Since  there  are  350  frequencies,  the  average  receiver  voltage 

-1/2 

modulation  depth  for  any  one  must  be  (350)  ~  0.042,  or  4.2 %  to  main¬ 

tain  a  constant  power -bandwidth  product.  This  assumes  a  100%  speckle 
modulation  depth,  which  will  occur  only  for  purely  diffuse  targets.  If  we 
take  a  50%  modulation  depth  since  realistic  targets  will  probably  have  some 
specular  returns,  the  average  modulation  depth  drops  to  2.1%.  This  value 
is  below  the  5.  6%  peak  receiver  modulation  produced  by  a  multidither  sys¬ 
tem  operating  with  a  ±20°  dither  modulation. 

A  further  consideration  that  we  have  not  taken  into  account  is  that  the 
spectral  components  of  speckle-induced  receiver  noise  will  have  randomly 
varying  phases  with  respect  to  the  dither  modulation.  The  servo  synchro¬ 
nous  detectors  will  thus  further  discriminate  against  this  type  of  noise.  We 
are  developing  an  analytical  demonstration  of  the  magnitude  of  this  discrimi¬ 
nation  on  Contract  F30602-76-C-0021 .  Some  of  the  details  of  the  early 
stages  of  this  analysis  are  given  in  Reg.  6  and  in  Section  IV-8  of  this  report. 


*A  more  realistic  speckle-induced  spectrum  might  be  a  uniform  triangular 
or  gaussian  distribution  with  a  50  kHz  bandwidth. 
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3. 


Conclusions 


These  brief  studies  seem  to  indicate  that  target- signature 
effects  are  not  going  to  be  as  severe  as  was  thought  earlier.  More  realistic 
targets  and  motions  need  to  be  studied,  however,  and  other  effects  such  as 
more  servo  channels  or  receiver  signal-to- shot  noise  need  to  be  investigated. 
The  presence  of  turbulence  or  blooming  may  also  have  some  effect. 
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IV. 


ANALYTICAL  TECHNICAL  ACCOMPLISHMENTS 


Much  of  the  analytical  work  performed  on  this  contract  has  been 
accomplished  using  computer  simulation.  The  simulation  codes  model  both 
the  details  of  a  COAT  servomechanism  (particularly  in  the  case  of  multi¬ 
dither  COAT)  and  the  details  of  the  propagation  of  a  coherent  beam  through  a 
turbulent  and  absorbing  medium.  The  codes  have  been  developed  to  their 
present  state  on  a  variety  of  programs  including  the  Hughes  IR&D  program 
and  several  contracts:  N60921 -74-C-0249  (NSWC),  F30602-73-C-0248 
(DARPA/RADC),  F30602-75-C-0001  (DARPA /RADC) .  The  results  of  the 
work  performed  on  this  contract  are  summarized  in  this  section. 

A.  Computer  Simulation  of  COAT  Compensation 

1  *  Summary  of  Capabilities  and  Limitations  of  the  Computer 

Simulation  Codes 

The  elements  of  the  adaptive  optics  computer  simulation  are 
indicated  in  the  flow  diagram  in  Fig.  27.  The  simulation  contains  computer 
routines  that  model  the  operation  of  an  adaptive  phase  compensation  system, 
the  propagation  of  an  optical  beam  through  a  turbulent  and  absorbing  medium, 
and  the  reflectivity  function  of  an  idealized  target.  The  adaptive  system 
simulator  consists  of  two  main  subblocks:  a  servo  routine  that  can  model 
either  phase  conjugate  (return-wave)  control  or  multidither  (outgoing-wave) 
control  and  a  phase  corrector  routine  that  can  model  segmented  or  continu¬ 
ously  deformable  mirrors  as  well  as  tracking  and  focus  control.  The  latter 
capability,  added  under  Amendment  No.  1  to  this  contract,  is  discussed  in 
detail  in  Ref.  6.  The  simulation  models  a  closed-loop  adaptive  system  in 
which  the  control  information  is  provided  by  the  detection  of  a  return-wave 
reflected  from  the  target. 

The  effect  of  the  atmosphere  on  the  laser  beam  propagation  is 
modeled  by  a  time-dependent  propagation  code  that  treats  the  effects  of  both 
turbulence  and  thermal  blooming.  The  time-dependent  code  is  used  on  both 
the  outgoing  and  return  paths.  The  difference  in  the  return  calculation  is 
that  the  propagation  is  modeled  by  a  low-power,  linear  code,  which  includes 
the  effects  of  the  inhomogeneities  caused  by  turbulence  and  the  absorptive 
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Fig.  27.  Adaptive  optical  system  simulation  flow  diagram.  Both 
multidither  and  phase  conjugate  control  algorithms  are 
simulated  as  well  as  atmospheric  turbulence  and  thermal 
blooming . 
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heating  introduced  by  the  high-power  transmitted  wave;  the  outgoing 
calculation  includes  nonlinear  heating  effects  (i.  e.  ,  the  outgoing  beam 
affects  the  medium). 

The  time-dependent  propagation  algorithm  used  in  the  simulation  is 
basically  a  multipulse  code  in  which  the  laser  energy  is  propagated  in  a 
series  of  short  pulses.  The  pulse  length  is  assumed  short  enough  that  no 
heating  occurs  during  the  pulse.  Furthermore,  it  is  assumed  that  the 
heating  caused  by  a  pulse  has  reached  a  steady  state  by  the  time  the  next 
pulse  arrives.  The  term  steady  state  is  used  here  in  the  sense  that  all 
density  gradients  traveling  at  the  sound  velocity  are  assumed  to  have  left 
the  region  occupied  by  the  laser  beam  by  the  time  the  next  pulse  arrives. 

The  only  gradients  present  are  those  that  are  convected  with  the  local 
crosswind  velocity.  Each  pulse  "sees"  the  medium  created  by  the  steady 
state  heating  of  all  previous  pulses  but  is  unaffected  by  its  own  heating. 

This  model  allows  us  to  simulate  either  pulsed  or  cw  propagation  depending 
on  the  time  between  pulses  relative  to  the  time  it  takes  an  air  parcel  to 
traverse  the  beam.  Continuous  wave  propagation  conditions  are  simulated 
by  allowing  significant  overlap  between  the  regions  heated  by  successive 
pulses.  On  this  contract,  we  have  studies  only  cw  propagation. 

The  use  of  a  time-dependent  propagation  code  differs  from  the 

approach  taken  by  Bradley  and  Herrmann  of  Lincoln  Laboratory,  who  have 

19 

used  a  cw  code  in  their  phase  compensation  studies.  Although  the  use  of 
a  time-dependent  propagation  code  greatly  increases  the  complexity  and 
cost  of  the  calculation,  we  have  felt  compelled  to  use  this  approach  for 
two  reasons.  First,  the  time-dependent  approach  is  closer  to  the  physics 
of  the  striation  that  we  are  trying  to  simulate  since  we  are  primarily 
interested  in  adaptive  optical  systems  that  respond  in  a  time  short  compared 
to  the  time  it  takes  the  atmosphere  to  change.  Such  adaptive  systems  per¬ 
form  an  essentially  instantaneous  determination  of  the  required  phase  front. 
Moreover,  the  determination  is  continuously  updated  as  the  medium 
changes,  rather  than  to  be  put  in  and  held  fixed  until  a  steady  state  is 
reached,  as  is  implicitly  assumed  when  a  cw  propagation  code  is  used. 

The  distinction  between  the  time-dependent  and  cw  approaches  is  espe¬ 
cially  important  for  the  simulation  of  adaptive  systems  that  employ 
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return- wave  algorithms.  Whereas  it  can  be  argued  that  a  cw  code  in 
conjunction  with  a  maximization  routine  closely  models  the  behavior  of  an 
adaptive  optical  system  employing  an  outgoing-wave  multidither  algorithm, 
it  does  not  appear  possible  to  model  return- wave  systems  with  the  cw 
approach  since  such  systems  do  not  explicitly  maximize  the  irradiance  on 
the  target.  Instead,  they  continuously  set  the  phase  of  the  transmitted  wave 
equal  to  the  conjugate  of  the  phase  of  the  field  returned  from  the  target. 
There  does  not  appear  to  be  any  reason  to  expect  that  the  final  state  reached 
by  such  a  system  can  be  obtained  from  a  cw  calculation  or  a  sequence  of 
such  calculations. 

The  other  reason  for  the  use  of  a  time-dependent  propagation  code 
in  our  studies  is  the  desire  to  follow  the  time  development  of  the  correction 
process  especially  in  the  case  of  adaptive  systems  employing  outgoing-wave 
multidither  control  algorithms.  An  outgoing- wave  multidither  control  sys¬ 
tem  attempts  to  maximize  the  irradiance  at  the  target  at  each  instant.  We 
wish  to  know  if  this  instantaneous  maximum  control  policy  leads  to  a  global 
maximum  in  the  thermal  blooming  correction  problem  or,  instead,  yields 
a  secondary  maximum.  This  can  be  determined  by  comparing  the  maximum 
obtained  with  the  time-dependent  and  the  cw  codes  because  the  cw  approach 
yields  the  global  maximum. 

The  time  sequence  of  events  in  the  simulation  is  shown  in  Fig.  28. 

A  high-power  pulse  of  laser  energy  is  propagated  to  the  target  and  the 
resulting  change  in  the  density  of  the  air  in  the  region  occupied  by  the  beam 
is  calculated,  taking  into  account  the  convective  effect  of  the  local  crosswind 
velocity.  It  is  assumed  that  the  adaptive  phase  compensation  is  applied  just 
prior  to  the  firing  of  the  next  high-power  laser  pulse.  Moreover,  we  are 
presently  assuming  that  the  adaptive  optical  system  has  infinite  bandwidth 
so  that  it  is  able  to  determine  and  apply  the  phase  correction  instantaneously 
without  any  delay  or  time-averaging  effects.  The  information  required  to 
determine  the  phase  correction  is  obtained  by  transmitting  a  low-power 
reference  pulse  to  the  target  and  reflecting  the  same  from  the  target.  This 
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Fig.  28.  Time  sequence  of  events  in  multipulse  adaptive 
optics  simulation: 

(1)  No  heating  during  a  pulse  (t  <<  t acous^c)> 

(2)  Steady  state  heating  for  t-t^  > 

(3)  Adaptive  optics  correction  applied  in  the 
interval  tdelay<  t-t.  <  tR. 
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reference  pulse  propagates  through  the  medium  that  will  be  seen  by  the 
next  high-power  pulse  and  thus  provides  the  necessary  phase  correction 
information.  This  process  is  then  repeated  for  each  subsequent  high-power 
pul  s  e . 

With  the  exception  of  some  runs  in  which  the  effect  of  the  atmosphere 
was  represented  by  a  single  phase  screen  in  front  of  the  transmitter,  all  of 
the  data  that  have  obtained  to  this  date  apply  to  the  simulation  of  adaptive 
systems  utilizing  a  return-wave  control  algorithm,  i.  e.  ,  systems  that  cor¬ 
rect  by  inserting  the  conjugate  of  the  phase  of  the  wave  returned  from  the 
target.  The  reason  for  this  concentration  on  the  phase  conjugate  simulation 
is  simply  that  such  simulations  are  significantly  easier  to  implement  and 
less  costly  to  run  than  outgoing- wave  simulations.  The  simulation  of  a 
return-wave  system  requires  simply  that  we  transmit  a  single  low-power 
reference  wave  to  the  target  and  a  single  low-power  reflected  wave  back  to 
the  transmitter  aperture.  The  phase  correction  is  then  equated  to  the  con¬ 
jugate  of  the  phase  of  the  reflected  field.  In  contrast,  an  accurate  multi¬ 
dither  simulation  requires  the  sending  and  receiving  of  many  low-power 
reference  pulses  for  each  high-power  pulse.  We  must  sample  the  modula¬ 
tions  of  the  integrated  irradiance  of  the  target  returns  at  a  rate  set  by  the 
highest  dither  frequency.  For  example,  it  was  found  in  our  earlier  multi¬ 
dither  servo  system  simulations  that  a  sampling  time  increment  of  5  \xsec 
is  required  for  an  18-channel  system  having  an  upper  dither  frequency  of 
32  kHz  and  an  overall  convergence  time  of  1  msec.  Hence,  approximately 
200  samples/convergence  time  were  required  in  those  simulations.  It  was 
practical  to  take  200  samples /convergence  time  in  the  earlier  work  because 
of  the  simplicity  of  the  propagation  model  that  was  used.  In  the  present 
case,  however,  the  propagation  model  is  quite  complicated  and  relatively 
costly  to  run  and  thus  it  is  clearly  not  practical  to  do  an  equivalent  simula¬ 
tion.  Additional  details  on  the  computer  simulation  codes  —  servo,  propaga¬ 
tion,  mirror,  and  target  —  are  given  in  Ref.  30. 

2.  Turbulence  and  Blooming  Compensation 

30 

We  have  demonstrated  elsewhere  that  the  effect  of  atmo¬ 
spheric  turbulence  on  the  ensemble-averaged  irradiance  distribution  can  be 
accounted  for  by  introducing  a  random  phase  screen  in  front  of  the  transmitting 
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aperture.  The  appropriate  phase  distribution  is  that  which  would  be 

produced  by  a  point  source  located  on  the  target/1'  This  implies  that  an 

ideal  adaptive  optics  system  should  be  able  to  completely  compensate  for 

the  effects  of  the  turbulence  by  introducing  the  conjugate  phase  distribution 

as  a  correction.^  The  degree  to  which  this  is  achieved  depends  on  whether 

the  target  has  a  point-glint  scatterer  that  provides  the  requisite  phase 

information  and  on  the  degree  to  which  the  adaptive  optics  system  can 

reproduce  the  desired  phase  distribution. 

30 

Previously  we  have  demonstrated  that  phase  conjugate  COAT  con¬ 
trol  with  a  deformable  mirror  can  correct  a  beam  focused  through  a  turbu¬ 
lent  atmosphere  back  to  within  70  to  95%  of  the  diffraction-limited  (free- 

30 

space)  target  irradiance.  We  have  also  demonstrated  analytically  that  a 

multidither  COAT  system  will  arrive  at  the  same  phase  correction  as 

a  phase  conjugate  system  when  correcting  for  linear  (not  power-dependent) 

distortions.  The  amount  of  correction  depends  on  the  range  (Z),  turbulence 

strength  (C^.),  wavelength  (\),  transmitter  diameter  (DT),  and  number  of 

deformable  corrector  mirror  actuators  elements  (N  ).  The  dependence  of 

a  6  30 

the  residual  phase  error  after  correction  is  described  by 

«j>^>  =  0.051  Z  D*/3  N‘5/6  (5) 

where  «b  >  is  the  residual  mean  square  phase  error  across  the  aperture 
e 

after  correction.  This  result  includes  an  approximation  for  the  shape  of  a 
deformable  mirror  surface.  The  effective  Strehl  ratio  is  given  by 


The  value  of  S  is  plotted  in  Fig.  29  versus  C^  for  several  cases  of  interest. 

*In  this  discussion,  we  assume  that  amplitude  scintillation  effects  are  small. 
If  they  are  not  small,  the  effect  of  atmospheric  turbulence  cannot  be  repre¬ 
sented  simply  by  a  phase  screen  but  rather  we  must  also  introduce  a  ran- 
dom  apodization. 

"^Subject  to  the  assumption  that  amplitude  effects  are  negligible. 
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-  STREHL  RATIO 


Na  =  NUMBER  OF  ACTUATORS 


Fig.  29.  Strehl  ratio  associated  with  atmospheric  turbulence  ver 
sus  number  of  mirror  actuators. 
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To  provide  a  check  on  the  above  analytical  work,  we  have  performed 
two  types  of  mirror  simulations.  In  each,  the  effect  of  a  turbulent  phase 
screen  placed  in  front  of  an  aperture  was  compensated  by  a  mirror  with  a 
finite  number  of  actuators.  The  average  Strehl  ratio  obtained  from  a 
sequence  of  five  independent  phase  screens  was  determined  as  a  function  of 
the  number  of  actuators.  In  one  of  the  simulations  the  deformable  mirror 
software  was  used  in  conjunction  with  a  sinusoidal  multidither  COAT  algo¬ 
rithm  of  the  type  used  previously  at  HRL  in  our  COAT  servo  system  studies. 
In  the  other  simulation,  the  mirror  was  modeled  by  a  segmented  mirror  with 
piston  and  tilt  control  on  each  segment.  The  piston  and  tilt  settings  for  each 
segment  were  determined  by  a  least  squares  fit  to  the  random  phase  surface 
over  the  segment. 

The  agreement  between  the  piston  and  tilt  results  and  the  theoretical 
predictions  obtained  from  eqs.  (5)  and  (6)  is  very  good  (see  Ref.  30),  which 
gives  us  confidence  in  the  theoretical  results  given  in  these  equations.  We 
believe  that  these  results  are  representative  of  those  that  would  be  obtained 
with  a  deformable  mirror  in  the  absence  of  the  2Ntt  problem  is  avoided. 

The  agreement  between  the  theory  and  the  deformable  mirror- 

multidither  COAT  simulation  results  is  reasonably  good  for  moderate  values 
2  2 

of  C  Z  but  is  poor  for  large  values  of  C  Z.  We  attribute  the  poor  results 
N  2  hi 

obtained  at  large  values  of  C^Z  to  a  2NiT-type  of  behavior.  A  similar  prob¬ 
lem  is  observed  with  the  focus  control  computer  simulation.  The  deform¬ 
able  mirror  simulation  that  was  used  in  these  runs  has  a  2Nit  correction 
loop  that  introduces  a  2tr  correction  whenever  the  phase  difference  between 
actuators  exceeds  four  radians.  The  intent  is  to  suppress  2tt  errors  intro¬ 
duced  by  the  servo  system.  However,  if  the  phase  distortion  that  is  to  be 
removed  by  the  mirror  changes  by  more  than  four  radians  between  actua¬ 
tors,  the  2tt  "correction"  introduced  by  the  2Nir  correction  loop  is,  in  fact, 
a  2 it  error  and  the  mirror  performance  is  corresponding  degraded.  This 
problem  could  be  avoided  by  removing  the  2Ntt  correction  loop  but  then  we 
would  be  faced  with  2Ntt  servo  errors.  One  way  to  avoid  this  problem  is 
to  use  more  actuators  so  that  the  phase  change  between  actuators  never 
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Fig.  30.  Strehl  ratio  due  to  Kolmogorov  turbulence 


exceeds  four  radians.  Another  way  would  be  to  design  the  mirror  so  that 
the  likelihood  of  2Nir  errors  is  reduced,  in  which  case  the  2Ntt  correction 
loop  could  be  eliminated.  More  work  clearly  needs  to  be  done  on  this 
problem. 

We  have  made  several  computer  runs  to  determine  turbulence  com¬ 
pensation  using  only  the  tracking  and  focus  controls.  Figure~30  shows  the 

observed  performance  with  tilt  control  alone.  The  relationship  between  the 

2 

quantity  p  in  Fig.  30  and  C-^Z  is  shown  in  Fig.  31  for  several  cases  of 
interest.  Figure  32  shows  a  summary  of  the  compensation  results  achieved 
with  BAT  controls  compared  with  those  produced  by  a  37-element  multi¬ 
dither  COAT  control  system.  The  deformable  mirror  used  in  all  the  simu¬ 
lation  runs  has  37  actuators.  The  assumed  propagation  scenario  is  also 

2  -111/3 

shown  in  the  figure.  For  reference,  C^Z  £  21  x  10  m  is  considered 

to  be  very  strong  turbulence.  Also  shown  in  the  figure  is  the  theoretical 

limit  given  by  eqs.  (5)  and  (6)  for  a  37-element,  deformable  mirror,  phase- 

conjugate  COAT  system  that  has  perfect  phase-sensing,  but  which  is  limited 

by  the  ability  of  the  mirror  to  provide  the  desired  phase  front. 

There  are  several  noteworthy  features  of  the  data  in  Fig.  32.  First, 

the  agreement  between  analysis  and  computer  simulation  for  the  37- element 

deformable  mirror  COAT  system  is  very  good.  The  slight  discrepancy  at 
2 

large  values  of  C^-Z  may  be  caused  by  improper  operation  of  the  simulation 
routine  designed  to  remove  2mr  ambiguities,  as  discussed  above.  Second, 
significant  correction  is  obtained  with  the  tracking  or  focus  controls  when 
used  separately,  but  the  two-axis  focus  control  provides  the  most  correc¬ 
tion.  In  fact,  when  the  tracking  controls  are  combined  with  the  two-axis 
focus  control,  little  increase  in  Strehl  ratio  is  produced  over  that  observed 
with  the  focus  control  alone.  We  attribute  the  lack  of  a  significant  tilt  com¬ 
ponent  in  the  simulation  (and  thus  little  correction  with  the  tracking  control) 
to  our  assumption  of  an  infinite  outer  scale  for  the  turbulence.  If  we  had 
used  a  finite  outer  scale,  we  expect  that  the  tracking  correction  would  have 
been  more  significant. 

The  significantly  worse  performance  of  a  spherical1'  focus  control 
2  -111/3 

above  C-^Z  =  40  x  10  m  is  also  interesting.  Finally,  not  unexpectedly, 

^Two-axis  focus  control  with  both  axes  driven  by  the  same  dither  and  cor¬ 
rection  signals. 
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AVERAGE  GLINT  IRRADIANCE,  %  OF  FREE  SPACE 
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♦  37-ELEMENT  MULTIDITHER  CONTROL 


Cn2Z,  I0~"  ml/3 


Turbulence  compensation  produced  by  computer  simulation 
using  37 -element  deformable  mirror.  Multidither  COAT 
control  is  compared  to  BAT  (tracking  and  focus  control, 
focus  control  only  (both  spherical  and  2-axis),  and  track¬ 
ing  control  only.  An  infinite  turbulence  outer  scale  is  used 
here,  in  contrast  to  Fig.  30. 


a  37-channel  multidither  COAT  system  produces  significantly  more 
correction  than  the  BAT  controls  alone;  the  BAT  controls  produce  roughly  a 
factor  of  2  improvement  in  Strehl  ratio  and  the  37-channel  COAT  system 
can  improve  the  Strehl  ratio  another  factor  of  2. 

g 

a.  Turbulence  and  Blooming  —  As  we  reported  earlier, 

the  effects  of  turbulence  and  thermal  blooming  can  be  considered  separately, 
and  then  combined  in  a  mean- square  sum  fashion  to  obtain  the  beam  Strehl 
ratio  when  both  effects  are  considered  simultaneously.  With  this  result,  it 
is  not  too  surprising  that  COAT  compensation  for  turbulence  is  unaffected 
by  the  presence  of  thermal  blooming.  The  experimental  data  presented  in 
Section  III  demonstrate  this  conclusion.  Our  computer  simulations  also  con¬ 
firm  that  good  turbulence  compensation  is  possible  even  if  little  or  no  bloom¬ 
ing  compensation  occurs.  Two  examples  of  the  simulation  results  are  shown 
in  Figs.  33  and  34.  The  transmitted  beam  in  this  case  is  a  gaussian, 
truncated  at  the  10%  intensity  radius. 

There  are  two  principal  observations  to  be  made  on  the  data  in 

Figs.  33  and  34.  First,  there  is  very  little  thermal  blooming  compensation. 

This  result  is  not  consistent  with  our  experimental  work  nor  with  other 

1 6  1 7  19 

experimental  ’  or  computer  simulation  work.  We  feel  that  the 

absence  of  correction  in  our  simulations  is  caused  by  (1)  the  use  of  a  return- 

wave  correction  algorithm  and  (2)  by  an  isoplanatic  problem  associated  with 

30 

finite-width  glints.  Our  experimental  data  were  taken  with  a  fast  outgoing- 
wave  multidither  system.  The  data  obtained  by  others^’ were  also 
produced  by  outgoing- wave  systems,  but  in  each  case,  the  system  was 
effectively  a  slow  dither  system. 

We  are  beginning  to  accumulate  evidence  that  the  inability  of  a  phase 
conjugate  system  to  correct  for  blooming  —  and  its  tendency  to  decrease  the 
target  irradiance  with  strong  blooming  —  may  be  common  to  all  "return- 
wave"  COAT  systems.*  A  return- wave  COAT  system  is  one  that  operates 
only  on  the  wavefront  returned  from  the  target  to  derive  the  servo 


-  27 

❖We  have  observed,  as  have  others,  that  for  transmitter  powers  exceed¬ 
ing  the  optimum  power  for  blooming,  a  phase  conjugate  COAT  correction 
is  worse  than  doing  nothing  at  all. 
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Fig.  33.  Phase  conjugate  COAT  correction  for  thermal  blooming  at 
10.6  [im  for  a  truncated  gaussian  beam  with  beam  slewing. 
Cases  with  and  without  turbulence  are  shown. 


3719-63R2 


Fig.  34.  Phase  conjugate  COAT  correction  for  thermal 

blooming  and  turbulence  at  3 . 8  \xm  for  a  truncated 
gaussian  beam  with  beam  slewing. 
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correction  signals.  Phas  e- conjugate  systems,  return-wave  multidither 
systems/1'  and  hybrid  MTRIM-COATM  systems**  are  all  return-wave  systems. 
In  fact,  of  the  presently  known  adaptive  systems,-  only  the  multidither, 
outgoing-wave  COAT  system  appears  to  be  a  true  maximization  system. 

Our  simulation  work  has  also  indicated  that  referencing  off  of 
extended  target  regions  in  the  presence  of  blooming  will  give  incorrect 
correction  results.  The  problem  is  the  thermal  blooming  equivalent  of  the 
finite  isoplanatic  patch  size  encountered  in  trying  to  image  extended  targets 
through  turbulence.  If  the  target  reference  region  exceeds  an  atmospheric 
coherence  diameter, t  a  COAT  system  will  not  be  able  to  obtain  the  correct 
phase  errors  for  the  focused  propagation  path.  We  have  observed  during 
our  computer  simulation  work  with  phase  conjugate  correction  for  blooming 
that  this  problem  occurs  even  with  glints  less  than  one-half  the  diameter  of 
the  diffraction-limited  transmitted  beam  on  the  target.  We  plan  to  study 
this  problem  further  on  other  related  programs. 

The  second  observation  on  the  data  in  Figs.  33  and  34  is  that  the 
ability  to  compensate  for  the  effects  of  turbulence  is  apparently  not  degraded 
by  the  presence  of  thermal  blooming  effects.  The  net  improvement  in  beam 
quality  with  phase  conjugate  compensation  is  most  impressive  for  the  3.8  \xm 
case  in  Fig.  34  and  for  cases  even  at  10.6  |Jim,  with  large  transmitter  aper¬ 
tures.  This  is  to  be  expected,  of  course,  since  turbulence  effects  are  most 
pronounced  in  the  cases  of  large  beams  and  shorter  wavelengths. 

B.  Complex  Targets  and  Speckle  Modulations 

As  part  of  the  contract  redirection  defined  by  Amendment  No.  2, 
we  have  been  pursuing  a  study  of  the  behavior  of  a  multidither  COAT  system 
when  it  is  subjected  to  multiplicative  noise  in  the  receiver.  The  type  of 


*The  reflected  coherent  light  sensed  by  the  receiver  is  dithered,  rather 
than  the  transmitted  beam. 

**  Transmitting  -Imaging:  broad  band,  incoherent  target  returns  are  sensed 
by  an  image-compensating  Imaging-COAT  system  and  this  information  is 
used  to  correct  the  phase  of  the  transmitted  beam. 

t The  distance  between  two  points  in  the  observation  plane  for  which  rays 
emanating  from  a  point  source  see  the  same  optical  path  length. 
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noise  considered  is  very  similar  to  the  modulations  produced  by  complex- 
target  speckle  effects  (see  Section  III  and  Ref.  29).  Our  analytical  studies 
have  utilized  both  closed  form  analysis  and  computer  simulation  of  an 
18-channel  multidither  servo. 

1 .  Analytical  Models  for  Speckle  and  COAT  Servo 

In  the  previous  quarterly  report,  ^  we  developed  a  statistical 
analysis  intended  to  predict  and  explain  the  behavior  of  a  COAT  system  sub¬ 
jected  to  multiplicative  speckle  modulations.  This  analysis  will  now  be 
used  to  calculate  expectation  values  of  various  COAT  system  parameters. 
We  then  compare  these  calculations  to  data  gathered  from  computer 
simulations. 

The  first  parameters  of  interest  are  P-^,  the  average  power  of  the 
dither  modulation  signal,  and  Pg,  the  average  power  of  the  speckle  signal 
within  ±400  Hz  of  any  dither  frequency.  The  expression  derived  in  Ref.  6 
for  P^  for  N  servo  channels  is 


n 


7 


where  K  is  a  receiver  proportionality  constant,  i[i  is  the  dither  amplitude, 
and 


N 

b„  s  12  Sin(|3n  ‘Pm1  '  (7) 

m=l 
n^m 

The  J  (40  and  J,  (41)  functions  in  eq.  (7)  are  the  zero-  and  first-order  Bessel 

O  1  ,|_j^ 

functions  and  is  the  phase  error  signal  in  the  i  control  channel.  The 
expectation  value  of  P^  is 
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k2  /4J0(t)  i2 

T 


<PD>  -  T-1-W„2X  1  E<b^ 


N 


n 


(8) 


n=l 


Since  t>n  is  a  function  of  the  random  variables  (3^,  ••  •>  (3j^>  it  is  also 

a  random  variable  whose  expectation  value  is  expressed  as 


dpNP(Pl)P(p2)..  p(pN)  b^pj,  p2---  pN> 


(9) 


where  P  (P^)  is  the  probability  density  of  p  . 

We  assume  that  the  p.'s  are  uniformly  distributed  between  the  phase 
angles  qq  and  aQ  +  a  as  indicated  in  Fig.  35.  Since  a is  arbitrary,  it  is 
convenient  to  define  it  as  a  aQ  =  a/Z’,  this  will  make  the  integrals  in  eq.  (9) 
easier  to  calculate.  Physically,  a  is  the  instantaneous  mean  phase  error, 
which  the  COAT  system  is  trying  to  drive  to  zero.  With  this  definition,  as 
a  approaches  zero,  the  COAT  array  converges.  Conversely,  as  a  approaches 
2r,  the  array  approaches  a  convergence  level  of  zero. 

For  any  given  (3. ,  the  probability  density  is  P(p^)  =  a~ *  over  a  region 
-a/Z  <  a.  <  a/Z,  and  zero  elsewhere.  Using  the  above  information  plus  the 
definition  of  bn  in  eq.  (7),  we  can  write  eq.  (9) 

N  N  ra/Z  /•<*/ Z  /•tt/2 

<bn>  =  <'"3  H  /  dPm/  dPK  /  dpn  Sin  ,(in  -  Pm>  sin  (Pn  -  Pk> 

m=l  k=l  J-a/Z  J-a/Z  J-a/Z 

m£k,  n 
k£n 

N  rot/  Z  r<*/Z 

+  “'2m  I  |  dP„  311,2  'Pn-Pm*  •  <10> 

m=l  J-a/Z  J-a/Z 

m£n 
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Fig.  35.  Phaser  diagram  of  a  partially  con¬ 
verged,  N-element  COAT  array 
defining  the  quantities  a  and  a 0  used 
in  the  text. 
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Carrying  out  the  integrals  and  summations  in  eq.  (10)  gives 


<b2>  =  i  (N2  -  3N  +  2) 
n  2 


.  2 
sin 


(a/2) 


®  b  - 


+  i<N- 1: 


(ID 


Substitution  of  eq.  (11)  into  eq.  (8)  gives  the  expectation  value  of  the  average 
dither  power: 


<PD>  =  K 


* 

•  2/o\ 

[joM  J^)]2 

,/l 

_  J_+ J_\ 

Sm  (2)  I 

1  sin  al 

llN 

N2  N  / 

(a/2)2  ' 

L1  ■  «  J 

/I  1  \  f ,  sin2  a 

U2"nV  l  ■  «2  J; 


(12) 


The  expectation  value  of  P  will  now  be  calculated.  As  derived  pre¬ 


viously, 


P  = 


N  M 

7  k2  Xlaj2k 


(13) 


k=l  j=l 

where  1^  is  the  mean  convergence  level  and  a^.k  is  a  random  amplitude 
variable  that  characterizes  the  amplitude  of  the  spurious  modulation  at  fre¬ 
quency  We  have  assumed  that  there  are  M  discrete  modulation  fre¬ 

quencies  around  each  of  the  N  dither  frequencies  so  that  the  spurious  modula¬ 
tion  has  the  form 
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N  M 


M 


=  i  +  2ZDa  jkhkt  +  V  ’ 

k=l  j=l 


(14) 


where  4>.k  is  a  randomly  chosen  phase  between  zero  and  2tt. 

J  7 

Since  Ij^  is  a  function  of  the  independent  random  variables  (3^, 

•  •  •  Pjyf  and  the  aj^  are  themselves  independent  random  variables,  the 
expectation  value  of  Pg  can  be  expressed  as 


<P  > 
s 


N  M 

1 *2  <4>E  E 


2 

<a..  > 
jk 


(15) 


k=l  j=l 


In  our  computer  simulations,  the  a^  were  chosen  randomly  between  the 
values  of  0  and  amax-  The  probability  distribution  is  this  uniform  over  that 
interval  and  zero  elsewhere.  The  expectation  value  of  their  square  is 


2 

<a.,  > 
jk 


max 


max 


2, 

a  da 


max 


(16) 


The  mean  convergence  level,  1^,  is  calculated  to  be 


t2  , . ,  N  N 

(40  v—' 


M 


n=  1  m=  1 
n^m 


(17) 


2  . 

The  expectation  value  of  is 


<IM>  ”  /  d(3 1  /dP 2 


fdPNP(Pl)  P(P2)  **•  P(PN)IM(Pr  P2* 


*N> 


(18) 


93 


Substituting  the  definition  of  I  given  in  eq.  (17)  and  the  probability  densities 
(P  (P^)  =  d  *)  into  eq.  (18),  we  can  solve  the  expectation  value  of  1^..  These 
steps  will  be  omitted,  however,  for  brevity.  The  result  of  the  calculation  is 


Figure  3  6  plots  the  root-mean- square  value  of  I  found  from  eq.  (19)  as  a 
fraction  of  a  for  an  18-channel  COAT  system.  When  eqs.  (16)  and  (19)  are 
substituted  into  eq.  (15),  the  expectation  value  of  the  average  power  in  the 
speckle  modulation  signal  is  found  to  be 


Expectation  values  of  many  other  parameters  can  be  calculated  in  a 
similar  manner.  Each  calculation  will  give  expressions  that  are  functions  of 


N,  ijj,  and  a  and  possibly  other  variables  such  as  a  ,  K,  etc.  .  which  are 

max 

included  in  the  definitions  of  a  given  parameter.  For  any  particular  COAT 
system,  most  of  these  variables  will  be  specified.  The  variable  a,  how¬ 
ever,  is  a  statistical  parameter  which  characterizes  the  convergence 
level  based  on  an  intuitive  model.  It  is  therefore  questionable  whether  one 
could  ever  measure  such  a  parameter  experimentally.  Nevertheless,  know¬ 
ing  several  expectation  values  as  functions  of  allows  us  to  derive  expres¬ 
sions  for  one  expectation  value  as  a  function  of  another:  <Pq>  as  a  function 
of  <IM>,  f°r  example.  This  kind  of  analytical  prediction  can  be  checked 
against  the  real  world  and  is  therefore  more  useful. 

Figure  37  shows  plots  of  <P  >  and  <P_.>  versus  <L  r>  calculated  in 

s  D  M 

this  manner.  Several  curves  of  <P  >  corresponding  to  different  values  of 

s 

amax  are  s^own*  One  can  appreciate  how  speckle  modulations  effect  a 

COAT  system  by  considering  one  value  of  a  ,  e.  g. ,  a  =  0.  006. 

max  max 

Below  a  convergence  level  of  approximately  0.  5  the  power  in  the  dither 
signal  is  greater  than  the  power  in  the  speckle  signal  while  the  converse  is 
true  above  this  level.  Based  on  the  simulation  data,  the  COAT  system 
apparently  has  no  trouble  converging  to  L^.  =  0.  5,  where  power  ratio, 
p  =  <Pj-j>/<Pg>,  is  equal  to  unity.  However,  when  a  convergence  level  of 
approximately  1.  77  is  reached,  the  convergence  process  stops  and  a  random 
oscillation  between  the  convergence  level  of  0.  5  and  0.77  begins  (the  com¬ 
puter  simulation  runs  discussed  below  exhibit  this  behavior). 

From  the  data  in  Fig.  37  it  is  evident  that  for  any  given  a  the 

max 

ratio  Rj,  of  the  dither  power  to  the  speckle  power  varies  over  a  large 
range  of  values  as  a  function  of  convergence  level.  For  the  example  of 
amax  =  °*  power  ratio  is  very  small  at  a  convergence  level  of 

O.  77  (the  actual  value  is  about  0.  3).  With  the  speckle  power  so  much 
larger  them  the  dither  power  it  is  reasonable  to  expect  convergence 
problems. 

2.  Computer  Simulation  Code 

Along  with  this  analysis,  we  have  been  making  computer 

runs  using  the  detailed  COAT  servo  simulation  shown  in  Fig.  38.  The 

spurious  modulation  function  M  (as  defined  by  eq.  (14))  in  our  work  to  date) 

s 
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Fig.  37.  Theoretical  dither  power  and  speckle  power  versus 
convergence  level  for  an  18-element  COAT  array 
( vjj  =  20°). 
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Fig.  38.  COAT  servo  simulation  showing  addition  of  "speckle11  noise  generator 


is  introduced  as  multiplying  the  normal  receiver  signal.  For  the  simulation 

studies  performed  to  date,  we  have  chosen  M  =  21  and  spaced  these  noise 

frequencies  ±83  Hz  around  each  dither  frequency.  This  is  a  rather  arbitrary 

choice,  but  it  does  ensure  that  all  spurious  modulation  is  well  within  the 

400  Hz  servo  bandwidth.  An  amplitude  limitation  was  also  included  in  the 

simulation  to  keep  the  return  signal  from  exceeding  an  absolute  value  of  unity 

just  before  entering  the  loop  gain  amplifier  (noted  as  in  Fig.  37).  This 

modification,  which  avoids  instabilities  resulting  from  the  automatic  gain 

control,  is  used  in  the  laboratory  model  of  the  COAT  array.  Simulations 

were  run  at  three  values  of  broad-band  signal-to-noise  ratio  (receiver  shot 

noise):  1010,  40,  and  20.  When  (s/N)ghot  =  1010  and  amax  =  0  (no  sPeckle 

noise),  the  simulation  produces  a  time  history  of  the  power  on  the  target 

glint  as  shown  in  Fig.  39(a).  When  the  speckle  noise  modulation  shown  in 

Fig.  39(c)  is  present  (a  =  0.  008);  the  time  history  in  Fig.  38(b)  results. 

max 

In  this  case,  the  rms  value  of  M  is  only  0.  064  (around  M  =  0.  50),  but  the 

s  s 

maximum  peak-to-peak  value  is  about  0.  4. 

The  behavior  shown  in  Fig.  39(b)  for  a  4  0  is  typical.  In  general, 

max 

the  speckle  noise  has  two  effects:  (1)  the  maximum  percent  convergence  level 
is  lowered,  and  (2)  a  random  oscillation  of  the  intensity  occurs  between  the 
new  maximum  convergence  level  and  a  minimum  level.  This  oscillation  is 
much  slower  than  the  dither  oscillations  and  its  magnitude  varies  with  con¬ 
vergence  level.  Figure  40  shows  a  plot  of  the  maxima  and  minima  recorded 
during  these  simulations  as  a  function  of  a  .  The  open  loop  gain,  G  ,  was 
the  same  for  each  case.  At  high  average  convergence  levels  (low  values  of 

a  ).  the  random  oscillations  have  much  smaller  magnitudes  than  at  the 
max 

low  convergence  levels  (high  values  of  a  ). 

XXl«*X 

3 .  Analysis  Compared  with  Computer  Simulation 

When  we  examine  other  speckle  power  curves  in  Fig.  36 

which  correspond  to  different  values  of  a  and  cross  check  these  with  the 
r  max 

computer  simulation  data,  it  becomes  apparent  that  the  power  ratio  is  a  very 
important  parameter,  defined  as 


P 


<PD> 


<P  > 
s 


(21) 
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Fig.  39.  COAT  computer  simulation  time  histories. 

(a)  Broadband  shot  noise  signal- to-noise 
ratio  =  (S/Nshot  =  10*0  and  amax  =  0  (no 
speckle  noise),  (b)  amax  =  0.008. 

(c)  Speckle  noise  modulation  used  to 
produce  (b). 
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Then  as  the  convergence  process  begins,  the  value  of  p  is  much  greater  than 
1.  As  the  convergence  level  increases,  the  value  of  p  gets  progressively 
smaller.  For  a  given  amax>  the  COAT  system  should  converge  to  the  level 
corresponding  to  p  =  1  (where  the  <Pg>  and  <-^d>  curves  Fig.  35  cross). 

The  final  convergence  level  corresponding  to  p  =  1  and  that  corresponding 
to  some  smaller  value  of  p. 

We  can  use  the  data  in  Fig.  40  along  with  the  analysis  developed 
above  to  verify  our  hypothesis  that  the  convergence  level  oscillates  between 
two  fixed  values  of  p  and  to  find  these  values.  This  is  done  by  first  choosing 

p  (e.  g. ,  p  =  1).  Equation  (12)  is  then  used  to  calculate  <Pn>  as  a  function 

2  1/2  U 

of  a  (or  equivalently,  <1  >  ).  Equations  (20)  and  (21)  are  then  used  to 

2  1/2 

calculate  a  as  a  function  of  <1  >  .  Using  this  method,  we  have  found 

max  m 

that  o  =1.0  and  p  .  =  0.  3  correspond  to  the  maximum  and  minimum 

Kmax  rmin  r 

convergence  levels.  Figure  41  compares  the  results  of  our  analytical 
calculations  with  the  computer-generated  data  of  Fig.  39.  The  agreement 
between  the  analytical  theory  and  the  simulation  data  is  excellent  for  signal- 
to-noise  (shot  noise)  ratios,  S/N,  between  1010  and  40.  The  agreement  is 
not  as  good  for  the  S/n  =  20  case,  but  this  is  not  unexpected  since  the  com¬ 
puter  simulations  have  always  shown  degradation  in  COAT  system  perfor¬ 
mance  when  S/N  was  chosen  less  than  40,  with  or  without  speckle  modulations. 
These  data  thus  appear  to  verify  the  original  hypothesis  that  motivated  much 
of  the  statistical  analysis,  namely  that  the  maximum  and  minimum  conver¬ 
gence  levels  correspond  to  constant  values  of  the  dither-to- speckle  power 
ratio. 

The  good  agreement  between  our  simplified  analysis  and  the  computer 
simulation  data  has  given  us  confidence  that  we  can  accurately  predict  the 
performance  of  a  multidither  COAT  system  in  the  presence  of  spurious 
modulations.  Our  work  to  date  has  been  with  a  particular  COAT  system 
and  a  hypothetical  model  for  the  spurious  modulations  (chosen,  however,  to 
present  maximum  difficulty  to  the  COAT  system).  We  expect  that  the  analy¬ 
sis  can  be  extended  in  a  straightforward  manner  so  that  it  is  independent  of 
the  speckle  modulation  model,  depending  only  on  the  statistics  of  the  modula¬ 
tion  spectrum.  A  contract  that  is  now  in  progress  (F30602-76-C-0021)  will 
accomplish  this  extension.  We  also  expect  that  modifications  to  the  servo. 
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Fig.  40.  COAT  convergence  levels  as  a  function  of  speckle  noise  ampli¬ 
tude  as  observed  by  computer  simulation. 
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which  are  not  accounted  for  in  the  analysis  (the  details  of  the  low-pass  filter, 

?  \  it 

for  example),  will  affect  the  values  of  p  and  p  .  and  the  curves  <1  > 
c  rmax  "min  m 

versus  a  for  a  given  p.  The  ability  of  the  analysis  to  predict  the  con- 
niax 

vergence  levels  as  a  function  of  a  (or  some  other  measure  of  the  speckle 

m  ax 

noise)  should  be  unaffected. 

At  this  point,  it  appears  that  the  speckle  modulation  problem  can  be 
viewed  as  a  signal-to-noise  (S/N)  problem.  In  addition,  the  changes  in  COAT 
system  performance  caused  by  speckle  "noise"  are  independent  of  broadband. 
The  signal-to-noise  (shot,  thermal,  etc.  )  ratio  as  long  as  this  latter  S/N 
ratio  is  sufficient  for  stable  operation  the  absence  of  speckle  modulations. 

The  COAT  system  signal-to- speckle  noise  ratio  varies  dramatically  with 
convergence  level,  however,  since  the  speckle  noise  within  the  ditherband 
competes  so  strongly  with  the  COAT-generated  convergence  signals.  With 
this  viewpoint,  any  changes  in  the  COAT  servo  that  act  to  improve  the  S/N 
will  reduce  the  deleterious  effects  of  speckle  modulations.  For  example, 
we  expect  that  increasing  the  dither  amplitude  or  the  amount  of  low-pass 
filtering  will  allow  operation  at  higher  speckle-noise  levels. 
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V. 


SUMMARY  AND  CONCLUSIONS 


A.  Turbulence  and  Blooming  Compensation 

This  program  has  continued  to  build  the  experimental  capabilities 
of  the  DARPA/RADC  COAT  experimental  model  and  to  increase  its  ver¬ 
satility  as  a  breadboard  for  studying  multidither  COAT  concepts  and  per- 
formance,  This  system  has  proved  to  be  the  single  most  valuable  tool 
available  to  date  for  such  studies,  rivaled  only  by  the  detailed  computer 
simulations  that  have  been  developed  along  with  it.  The  COAT  system  now 
has  18  phase  control  channels,  two  tracking  control  channels,  and  one 
focus  control  channel.  For  propagation  studies,  there  are  several  thermal 
blooming  and  turbulence  generators,  including  a  100  m  outdoor  path  at 
Hughes  Research  Laboratories. 

The  DARPA/RADC  COAT  system  has  demonstrated  that  an  improve¬ 
ment  factor  of  about  1.5  is  possible  against  blooming  distortions,  with  the 
compensation  occurring  only  for  distortions  in  the  first  3  0%  of  the  propaga¬ 
tion  path.  Compensation  for  turbulence  and  fixed  distortions  is  not 
impaired  by  the  presence  of  thermal  blooming,  even  though  little  correction 
may  occur  for  the  blooming  itself.  This  observation  has  been  verified  by 

both  experiment  and  computer  simulation.  These  results,  along  with  those 

1-5 

from  an  earlier  program,  indicate  that  the  achievement  of  the  high 
focal-plane  irradiance  should  not  be  limited  by  atmospheric  propagation, 
provided  that  thermal  blooming  distortions  are  minimized  and  that  a  suit¬ 
able  COAT  system  is  used. 

Computer  simulation  data  indicate  that  a  return- wave  COAT  system 
may  not  be  a  suitable  concept  when  thermal  blooming  is  present.  Such 
systems  can  remove  turbulence  and  fixed  distortions  just  as  a  multidither 
system  does,  but  may  become  unstable  and  reduce  the  focal  plane 
irradiance  if  the  blooming  becomes  too  severe. 

In  all  of  the  experiments  performed  to  date,  the  multidither  tracking 
and  focus  (BAT)  controls  have  not  increased  the  average  peak  beam  irradi¬ 
ance  significantly  over  that  produced  by  the  phase  controls  alone.  The 
only  exception  to  this  observation  occurred  with  artificial  turbulence 
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distortions  introduced  with  a  phase  screen.  With  a  different  artificial 
phase  screen,  the  BAT  controls  had  no  effect.  This  somewhat  unexpected 
observation  may  be  a  result  of  the  characteristics  of  our  artificial  tur¬ 
bulence  screen  since  real  turbulence  should  have  significant  tilt  compo¬ 
nents  (see  the  results  in  Section  IV  as  well  as  Ref.  11). 


Bo  Complex  Target  Effects  Caused  by  Speckle  Modulations 

Problems  with  COAT  system  convergence  caused  by  speckle- 
modulations  have  not  been  observed  with  an  extended,  scotchlite  target 
used  with  the  18-c'hannel  COAT  system.  Artificially-generated  amplitude 
modulations  have  affected  the  COAT  convergence  level,  however,  in  both 
experiments  and  computer  simulations.  It  remains  to  be  determined  whether 
real  targets  can  produce  significant  speckle  modulations  that  will  affect 
real  multidither  COAT  systems.  An  analytical  model  has  been  developed 
that  can  accurately  predict  the  convergence  level  and  stability  of  a  multi¬ 
dither  COAT  system  once  the  system  is  defined  and  once  the  statistics 
(such  as  the  power  spectral  density)  of  the  spurious  modulations  are  known. 
Additional  work  is  continuing  on  this  subject,  supported  by  DARPA/RADC 
Contract  F30602-76-C-0021. 


C.  Recommendations  for  Additional  Investigations 

The  work  on  this  contract  has  satisfied  the  original  goals  of  the 
program  and  has  contributed  substantially  to  the  understanding  of  COAT 
systems,  their  capabilities,  and  their  limitations.  There  are  several 
areas,  however,  where  additional  work  is  needed.  First,  we  need  to  deter¬ 
mine  what  targets  and  scenarios  can  produce  speckle  modulations  that  are 
strong  enough  to  degrade  multidither  COAT  performance.  This  work  will 
involve  analytical,  computer  simulation,  and  experimental  efforts  and 
should  include  some  experimental  work  at  10.  6  pm,  or  another  convenient 
infrared  wavelength.  A  current  contract,  F30602-76-C-0021,  is  aimed 
at  some  of  this  work,  but  does  not  include  10.6  pm  experimental 
measurements . 
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A  second  area  of  investigation  is  in  thermal  blooming  compensation. 
All  of  the  work  to  date,  at  Hughes  as  well  as  at  other  laboratories,  indi¬ 
cates  that  present  COAT  techniques  will  not  do  much  for  thermal  blooming 
with  cw  lasers.  Other  techniques  such  as  different  multidither  concepts, 
irradiance  tailoring,  wavelength  selection,  and  pulsed  operation  should  be 
pursued  to  find  those  concepts  and  systems  that  minimize  thermal  blooming 
or  increase  the  ability  of  COAT  to  compensate  for  those  blooming  dis¬ 
tortions  that  do  occur. 

There  is  still  much  to  be  done  in  the  area  of  computer  simulation. 

The  codes  developed  to  date  for  both  COAT  and  propagation  have  been  invalu¬ 
able  in  identifying  and  studying  propagation  phenomena  and  compensation 
for  them.  They  have  also  been  essential  for  designing  and  evaluating  COAT 
hardware.  Additional  work  is  required  to  improve  these  codes  as  well  as 
to  utilize  them  in  studying  COAT  compensation  with  various  system  and 
target  scenarios.  We  need  to  develop  efficient  codes  that  can  simulate 
thermal  blooming  compensation  with  multidither  COAT  system  and  realistic, 
multiple-glint  targets. 

Since  most  of  the  high-energy  laser  applications  of  COAT  are  at 
infrared  wavelengths,  much  experimental  work  needs  to  be  done  to  demon¬ 
strate  the  required  hardware  technology,  which  can  be  quite  different  with 
high-power,  infrared  lasers  than  with  low-power,  visible  wavelength  devices. 
The  HICLAS  program  (contract  N60921-76-C-0008  and  its  follow-on)  will 
address  some  of  these  technology  questions.  In  addition,  low-power  COAT 
experiments  with  infrared  wavelengths  are  in  order  to  investigate  target 
effects  such  as  speckle  and  glint  destruction  that  may  be  quite  different  at 
infrared  wavelengths  than  in  the  visible. 

Finally,  as  promising  new  COAT  servo  and/or  system  concepts  are 
invented,  they  should  be  investigated  in  low-power,  breadboard  experiments 
to  establish  what  advantages  they  may  have  in  compensation  performance, 
versatility,  simplicity  in  hardware  design  and  implementation,  or  relia¬ 
bility  of  operation  under  a  variety  of  operational  conditions. 
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